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ABSTRACT 

This  report  investigates  the  suitability  of  a  number  of  broadcast  band  trans¬ 
mitters  of  opportunity  for  a  passive  radar.  Specifically,  this  report  reviews  the 
design  of  active  radar  waveforms,  measures  emissions  for  representative  trans¬ 
mitters  of  opportunity,  and  assesses  the  suitability  of  these  signals  for  various 
tjrpes  of  radar  observation. 

This  report  shows  that  it  is  indeed  possible  to  construct  the  suggested 
passive  radar,  although  performance  is  limited  and  depends  largely  on  other 
parameters,  such  as  receiver  bandwidth,  antenna  pattern  and  observation  time. 
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Waveform  Analysis  of  Transmitters  of  Opportunity  for 

Passive  Radar 


EXECUTIVE  SUMMARY 

Most  radar  systems  work  by  actively  illuminating  the  target  environment  with  a  specific 
waveform  of  appropriate  energy  and  measuring  the  energy  scattered  by  the  target  to 
determine  target  parameters,  such  as  range,  bearing  and  velocity.  Accordingly,  most 
radar  systems  use  a  radar  waveform  that  optimises  performance  and  much  research  has 
been  completed  in  the  field  of  active  radar  waveform  design. 

An  alternative  is  not  to  deploy  a  dedicated  radar  transmitter,  but  rather  to  make  use  of 
transmissions  of  opportunity  such  as  radio  and  television  broadcast  signals.  Transmitter¬ 
less  passive  radar  has  a  number  of  advantages,  such  as  being  able  to  operate  without 
being  detected,  which  makes  it  more  attractive  than  conventional  active  radar  for  certain 
defence-related  missions.  It  also  has  disadvantages,  one  which  is  that  the  illumination 
waveform  is  almost  certainly  not  optimal  for  radar  operation. 

This  report  investigates  the  waveforms  generated  by  a  number  of  representative  trans¬ 
mitters  of  opportunity,  and  shows  how  these  waveforms  affect  the  performance  of  a  passive 
radar  system  utilising  them.  It  attacks  the  problem  of  opportunistic  passive  radar  with 
waveform  design  theory  and  an  assessment  of  the  suitability  of  the  transmissions  of  op- 
portimity  to  the  matched  filter  receiver.  Previous  analysis  of  this  kind  has  been  for  very 
specific  radar  systems  and  little  work  is  available  which  addresses  the  issues  with  radar 
waveform  theory. 

This  report  details  the  requirements  for  constructing  a  passive  radar  and  shows  how  its 
performance  varies  as  a  function  of  receiver  bandwidth,  antenna  pattern  and  observation 
time.  The  results  that  this  report  formulates  and  the  conclusions  it  draws  provide  a  solid 
foundation  for  anyone  desiring  to  construct  a  passive  radar. 
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Chapter  1 

Introduction 


Most  radar  systems  work  by  transmitting  a  specific  electro-magnetic  signal  into  a  region 
and  detecting  the  wave  scattered  by  an  object  of  interest.  By  analysing  this  scattered 
signal,  it  is  possible  to  determine  information  about  the  object,  such  as  range,  bearing, 
cross-section  and  velocity.  This  signal  is  referred  to  as  the  radar  waveform,  and  the  ability 
of  the  radar  to  perform  these  previously  mentioned  tasks  rests  heavily  upon  it. 

Accordingly,  most  radar  systems  use  a  radar  waveform  that  optimises  performance  and 
much  research  has  been  completed  in  the  field  of  radar  waveform  design  [1-4]. 

An  alternative  radar  system  would  be  one  which  uses  electro-magnetic  signals  that 
are  already  present  in  the  target  environment.  This  configuration  is  often  referred  to  as 
passive  radar,  as  the  transmitter  is  not  part  of  the  system. 

Requiring  no  transmitter,  passive  radar  has  many  benefits,  including  being  typically 
smaller  and  more  portable,  less  expensive,  requiring  less  design  efibrt  and  no  shielding  or  li¬ 
censing  issues  (for  spectrum  allocation).  Passive  radar  could  also  provide  great  advantages 
in  defence  as  the  operation  of  typical  active  radar  is  easily  detected  (the  transmitted  wave¬ 
form  is  very  distinct  and  usually  occupies  known  or  easily  determined  carrier  frequencies). 
For  this  reason,  active  radar  has  the  following  disadvantages: 

•  Ships,  submarines  and  planes  using  active  radar  for  navigation  or  target  detection 
may  find  it  impossible  to  travel  inconspicuously  within  particular  regions. 

•  It  may  prove  impossible  to  use  an  active  radar  to  inconspicuously  monitor  a  sensitive 
region  of  interest. 

•  Many  offensive  measures,  such  as  High-speed  Anti-Radiation  Missiles  (HARMs), 
work  by  following  the  enemy’s  radar  transmissions  to  their  soiurce  [5]. 

Transmitter-less  passive  radar  is  undetectable  and  thus  presents  none  of  the  above 
limitations. 

The  concept  of  utilising  transmitters  of  opportunity  is  not  new,  although  the  major¬ 
ity  of  available  literature  considers  only  the  possibility  of  using  illumination  provided  by 
another  radar.  Schoenenberger  and  Forrest  [6]  and  Green  [7]  both  detected  signals  from 
the  radar  transmitter  located  at  an  airport  to  determine  the  location  of  planes  arriving 
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at  and  departing  from  that  airport.  More  general  investigations  of  passive  radar  which 
utilise  transmissions  from  other  radar  have  been  perforined  by  Thompson  [8],  Munich 
and  Schecker  [9]  and  Lindop  [10].  The  advantage  of  such  radar  systems  is  that  the  radar 
waveform  is  optimised  for  target  detection  and  parameterisation. 

It  is  conceivable  that  a  passive  radar  employ  the  use  of  transmissions  from  more  general 
sources  such  as  radio  and  television  broadcast  towers.  In  this  system,  the  radar  waveform 
would  not  be  optimum,  but  an  effectively  random  one  generated  by  modulating  infor¬ 
mation  such  as  speech,  music  and  video.  Joint  research  between  the  University  College, 
London,  and  Britain’s  Defence  Research  Agency  [11—13]  showed  that  such  a  radar  was 
possible  using  television  broadcast  transmissions,  although  performance  was  very  limited. 
Also,  Sahr  and  Lind  [14]  demonstrated  a  radar  which  utihsed  transmissions  from  a  local 
FM  radio  station  to  monitor  the  state  of  the  upper  atmosphere.  Due  to  its  application  to 
defence,  it  is  believed  that  much  more  research  has  been  performed  for  this  class  of  passive 
radar  than  is  available  in  the  open  literature. 

This  report  investigates  the  waveforms  generated  by  such  transmitters  of  opportunity, 
and  shows  how  these  waveforms  effect  the  performance  of  a  passive  radar  system  utilising 
them.  The  work  presented  in  this  report  is  a  generalisation  of  the  work  performed  by 
the  authors  mentioned  in  the  previous  paragraph,  who  have  attempted  to  analyse  very 
specific  passive  radar  systems.  This  report  attacks  the  problem  of  opportunistic  passive 
radar  with  waveform  design  theory  and  an  assessment  of  the  suitability  of  the  transmissions 
of  opportimity  to  the  matched  filter  receiver. 

This  report  will  show  that  it  is  indeed  possible  to  construct  the  suggested  passive  radar, 
although  performance  is  limited  and  depends  largely  on  other  parameters,  such  as  receiver 
bandwidth,  antenna  pattern  and  observation  time. 

The  rest  of  this  chapter  completes  the  introduction  to  the  problem.  It  outlines  the 
contents  of  this  report  and  discusses  the  approximations  and  assumptions  used  throughout 
it. 

Chapter  2  presents  a  brief  theory  of  radar  waveforms.  The  importance  of  the  waveform 
to  radar  performance  has  already  been  stressed  and  this  chapter  explains  this  relation¬ 
ship.  It  introduces  the  matched  filter  and  ambiguity  function  and  gives  the  reader  some 
appreciation  of  radar  waveform  design. 

Chapter  3  provides  some  details  on  what  transmissions  of  opportunity  exist  and  pro¬ 
poses  a  radar  system  which  utilises  these  transmissions.  It  provides  a  picture  of  how 
a  non-cooperative  passive  radar  could  be  constructed  however  it  is  left  to  the  following 
chapter  to  determine  whether  it  would  work. 

Chapter  4  applies  the  theory  of  chapter  2  to  the  proposed  system  of  chapter  3.  It  shows 
that  the  suggested  radar  could  work  and  details  how  the  other  variables  of  the  radar  can 
be  adjusted  to  provide  the  best  possible  performance. 

Chapter  5  provides  a  number  of  simulated  examples  of  the  ideas  presented  in  this 
report.  Although  many  parameters  of  the  radar  and  environment  are  simulated,  actual 
waveforms  from  transmitters  of  opportunity  are  used.  These  simulations  provide  some 
confirmation  of  the  results  attained  in  chapter  4  and  an  initial  picture  of  the  performance 
of  the  suggested  passive  radar. 


2 


DSTO-TR-0809 


The  topic  of  waveform  analysis  for  passive  radar  is  a  vast  subject  and  this  report 
attempts  only  to  be  one  part  of  a  more  detailed  investigation.  Topics  of  study  which  build 
on  the  results  presented  in  this  report  and  suggestions  on  physical  implementation  of  the 
proposed  radar  system  are  given  in  the  concluding  chapter  of  this  report. 
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Chapter  2 

Radair  Waveform  Theory 


The  ability  of  a  radar  to  detect  targets  amongst  returns  from  other  objects  (clutter)  in  a 
region  of  interest  and  the  ability  to  determine  parameters  of  these  targets,  such  as  range, 
bearing,  size  and  velocity,  depends  largely  on  the  radar  waveform.  This  chapter,  and  the 
rest  of  this  report,  focuses  upon  the  effect  of  the  waveform  on  the  ability  to  detect  targets 
amongst  clutter  and  determine  the  range  and  velocity  of  these  targets.  The  reason  for  this 
is  that  the  accuracy  with  which  the  range  and  velocity  may  be  determined  is  controlled 
by  the  waveform  used  by  a  bistatic  radar  with  stationary  transmitter  and  receiver.  The 
ideas  presented  however,  can  be  simply  extended  to  extract  other  information  about  the 
target  such  as  acceleration  and  cross-range,  which  is  possible  in  synthetic  aperture  radar 
(SAR). 

Provided  the  wave  scattered  from  the  target  can  be  detected,  range  is  proportional  to 
the  delay  in  time  between  transmission  of  the  waveform  and  reception  of  the  scattered 
component,  and  velocity  is  proportional  to  the  shift  in  frequency  of  the  waveform  as  it  is 
scattered  from  the  target  (the  Doppler  shift).  Note  that  both  here  and  throughout  the 
report,  range  refers  to  the  total  distance  from  the  transmitter  to  the  target  to  the  receiver 
and  velocity  refers  to  the  rate  of  change  of  this  total  range.  Depending  on  the  location  of 
the  transmitter  and  the  receiver  and  the  direction  the  target  is  travelling,  these  parameters 
may  not  reflect  the  actual  range  and  velocity  of  the  target,  however  it  is  usually  possible 
to  calculate  actual  range  and  velocity  from  the  measured  values.  Chapter  3  details  the 
configuration  of  the  radar  being  investigated  in  this  report  and  shows  how  this  conversion 
would  be  performed  in  this  case. 

It  should  also  be  noted  that  throughout  this  report  the  Doppler  effect  is  approximated 
by  a  translation  in  frequency  with  no  distortion  to  the  signal’s  envelope,  that  is,  S'{f)  = 
S{f  +  a).  The  effect  on  a  signal  reflected  from  an  object  moving  at  a  constant  speed 
however,  is  actually  a  scale  in  frequency  and  the  corresponding  envelope  scale  necessary 
to  conserve  energy,  that  is,  S'{f)  =  {l//3)S{^f).  It  can  be  shown  [1]  that  the  above 
assumption  is  valid  if  {2v/c)Ta  <  1,  where  v  is  the  velocity  of  the  target,  T  is  thexadar 
“look”  or  integration  time,  a  is  the  amount  of  frequency  shift  and  c  is  the  velocity  of  light. 
This  criteria  is  readily  satisfled  throughout  this  report  for  the  expected  values  of  v,  T  and 
a. 


The  waveforms  used  in  early  radars  were  single  frequency  sinusoids.  These  were  trans¬ 
mitted  either  in  short  pulses  in  order  to  accurately  measure  return  time  delay  or  for  very 
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long  durations  in  order  to  accurately  measure  Doppler  shift.  Detection  of  the  scattered 
wave  was  ensured  to  a  particular  range  by  adjusting  the  power  in  the  transmitted  wave¬ 
form.  This  presented  a  problem  when  the  desire  was  to  accmrately  measure  range,  as  the 
short  pulse  duration  of  the  waveform  required  the  transmitters  to  cope  with  very  high 
peak  voltages  to  attain  the  required  power. 

In  1953  Woodward  [15]  considered  using  waveforms  of  more  than  a  single  frequency 
and  showed  that  the  ability  to  detect  target  range  and  velocity  was  actually  a  function  of 
waveform  shape,  not  pulse  duration.  He  laid  the  foundations  for  a  unified  theory  which 
relates  the  waveform  to  the  ability  of  the  radar  both  to  detect  targets  and  to  resolve  target 
range  and  velocity.  It  is  an  overview  of  this  theory  which  is  presented  in  this  chapter. 

Usually,  this  theory  is  used  in  waveform  design,  that  is,  determining  what  waveform  a 
radar  would  use  given  a  specific  desired  performance.  This  project  proposes  to  investigate 
waveforms  emitted  from  transmitters  of  opportunity,  so  this  ability  is  not  directly  impor¬ 
tant.  However,  understanding  what  is  necessary  in  designing  a  satisfactory  waveform  is 
important  in  determining  whether  a  given  waveform  will  perform  satisfactorily. 


2.1  Matched  filter  receivers 

Independent  of  which  waveform  is  used,  it  is  an  obvious  requirement  to  extract  the  most 
information  about  the  target  as  possible  from  the  received  signal.  This  involves  minimising 
the  inaccuracies  due  to  the  various  noises  which  invariably  corrupt  the  received  signal. 

The  receiver  which  is  used  by  the  vast  majority  of  radars  and  communication  receivers 
is  the  one  which  optimally  detects  the  transmitted  signal  in  the  presence  of  additive  white 
Gaussian  noise.  It  is  referred  to  as  the  matched  filter  and  takes  the  form: 

h{t)=s{T-t)  (1) 

Where  h{t)  is  the  impulse  response  of  the  matched  filter,  s{t)  is  the  signal  desired  to  be 
detected  and  T  is  the  delay  of  the  filter. 

This  result  is  derived  by  maximising  the  ratio  of  the  probability  that  the  received 
signal  contains  s{t)  to  the  probability  that  it  does  not,  a  ratio  referred  to  as  the  likelihood 
function.  The  same  result  can  be  determined  by  maximising  the  signal-to-noise  ratio 
(SNR)  of  the  output  of  the  filter.  This  is  the  method  by  which  North  [16]  originally 
derived  the  above  expression.  For  a  detailed  derivation  of  the  matched  filter  the  reader  is 
referred  to  either  Cook  and  Bernfeld  [1]  (pages  18-33)  or  Van  Trees  [17]  (pages  275-279). 

In  the  case  of  radar,  the  signal  to  be  detected  is  the  component  of  the  radar  waveform 
scattered  by  the  target.  As  it  is  desired  that  the  filter  optimally  detect  the  signal  scattered 
from  a  target  whose  location  causes  a  delay  in  time  r  and  whose  velocity  causes  a  Doppler 
shift  (f)  of  the  transmitted  waveform  Xo{t),  then: 

s{t)  =  Xo{t  +  (2) 

which  implies  that  the  optimum  filter  is: 

h{t)=Xo{T-t  +  r)e^'^^'^-^'^  (3) 
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Note  that  it  is  assumed  the  bandwidth  of  the  waveform  is  small  compared  to  its  carrier 
frequency.  It  is  for  this  reason  that  a  frequency  shift  of  Xo{t)  can  be  written  as  a  phase 
shift. 

The  output  of  the  matched  filter  at  time  T  is: 

M(T)  =  /  Xr{v)h*{T  -v)dv 

Jo 

=  r  Xr{iy)xl{T  -{T-u)  + 

Jo 

=  f  Xr{v)xl{u  +  du  (4) 

Jo 

where  Xr{t)  is  the  received  signal. 

The  delay  required  by  the  matched  filter,  T,  is  the  radar  integration  time  and  will  be 
assumed  a  constant  from  this  point  on.  Instead,  the  output  of  the  matched  filter  will  be 
expressed  as  a  function  of  t  and  (j),  which  corresponds  to  the  range  and  velocity  of  the 
target  which  is  being  tested  for  detection.  Adding  a  normalising  factor,  the  final  expression 
for  the  matched  filter  receiver  for  radar  is: 

M(T,4>)  =  ^j^Xr{i)x*oit  +  t)e'‘^^*dt  (5) 

In  practice,  the  matched  filter  is  often  implemented  on  a  digital  computer  and  Xr{t)  and 
Xo{t)  are  converted  into  the  digital  signals  Xr[n]  and  XoN  respectively,  where  x[n]  =  x{nTs) 
and  Ts  is  the  sample  period.  In  this  case,  the  equation  for  the  matched  filter  receiver  is: 

=  Xr[n]x*[n  +  r]e“*^  (6) 

n=l 

where  N  is  the  number  of  samples  in  the  integration  time  {TgN  =  T).  Issues  arising 
from  digitally  performing  matched  filtering  are  discussed  in  section  2.5. 

The  expression  for  the  matched  filter  receiver  can  be  viewed  as  a  correlation  between 
the  transmitted  signal,  delayed  in  time  by  r  and  shifted  in  Doppler  by  and  the  received 
signal.  For  this  reason,  the  matched  filter  is  often  referred  to  as  the  correlation  receiver. 

The  result  of  the  matched  filter  receiver,  M(r,  ^),  is  a  complex  number.  As  it  is  the 
magnitude  of  this  quantity  which  provides  the  likelihood  of  a  target  being  present,  the 
following  strategy  for  target  detection  is  usually  adopted: 

•  Using  the  received  signal,  Xr{t),  a  copy  of  the  transmitted  signal,  Xo{t),  and  equa¬ 
tion  (5),  calculate 

•  Either  digitally  or  using  a  square-law  device,  calculate  |M(t,</i)P. 
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•  Compare  this  quantity  to  a  single  threshold,  1.  A  value  of  |M(t,  (^)P  greater  than 
I  indicates  the  presence  of  a  target  with  a  location  and  velocity  corresponding  to  a 
time  delay  r  and  Doppler  shift 

Van  Trees  [17]  showed  that  if  the  noise  corrupting  the  received  signal  is  additive  white 
Gaussian,  the  above  strategy  provides  optimum  detection. 

Another  point  of  note  is  that  equation  (5)  resembles  the  Fourier  transform  of  the 
product  of  the  signals  Xr{t)  and  x*{t  +  t)  and  equation  (6)  resembles  the  discrete  Fourier 
transform  (DFT)  of  the  product  of  the  signals  Xf\n\  and  x*[t  +  rj.  This  is  unsurprising 
as  correlating  against  different  frequency  shifts  is  merely  convolution  in  the  frequency 
domain.  This  operation  can  be  implemented  as  a  product  of  the  two  signals  in  the  time 
domain  and  the  corresponding  Fourier  transform. 

A  consequence  of  this  is  that  equation  (6)  could  be  implemented  using  a  fast  Fourier 
transform  (FFT)  algorithm.  In  practice,  this  and  other  techniques,  such  as  incorporating 
information  about  the  waveform  into  the  above  equations,  are  usually  required  for  the 
matched  filter  to  operate  in  real-time. 

Appendix  B.l  contains  Matlab  code  which  implements  the  matched  filter  receiver. 


2.2  The  ambiguity  function 


Let  x{t)  be  the  radar  waveform.  If  there  exists  a  single  target  in  the  region  of  interest, 
the  received  signal  is  ideally: 

Xr{t)  =  ax{t  -  (7) 

where  the  time  delay  td  is  due  to  the  location  of  the  target,  the  Doppler  shift  Wd  is  due 
to  its  velocity  and  a  represents  the  attenuation  of  the  transmitted  signal  due  to  path  loss 
and  scattering.  Prom  equation  (5),  the  result  of  the  matched  filter  receiver  and  square- law 
device  is: 


where: 


1  7^  . 

—  ax{t-td)e 

J—td 


a 


^  r  ^%{t')x*{t'  +  {t  +  td))e-^^^+^^'>^' dt' 

J-td 


\x{r  +  td,(f>  +  Wd)\^ 


(8) 


x(t,  ^  x{t)x*{t  -f  T)e 


(9) 
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The  function  x('r,  <p)  is  the  ambiguity  function.  It  describes  how  correlated  a  signal  is 
with  itself  delayed  in  time  by  an  amount  t  and  shifted  in  Doppler  an  amount  <f>.  For  this 
reason  it  is  also  referred  to  as  the  time-frequency  autocorrelation  function. 

When  the  signal  being  examined  is  discrete,  the  ambiguity  function  takes  the  form: 

x(t,  (f>)  =  ^Y,  (10) 


Appendix  B.2  contains  Matlab  code  which  implements  this  equation. 

It  should  be  noted  that  the  exact  definition  of  the  ambiguity  function  varies  throughout 
the  literature.  Often,  for  example,  the  term  ambiguity  function  is  used  to  refer  to  the 
quantity  |x('r,<^)|^. 

Equation  (8)  gives  some  initial  insight  to  the  significance  of  the  ambiguity  function: 
that  the  result  of  the  matched  filter  receiver,  when  the  signal  received  is  the  scattering  from 
a  single  target,  is  the  ambiguity  function  of  the  transmitted  signal,  scaled  and  shifted  to 
be  centred  on  the  time  delay  and  Doppler  shift  corresponding  to  the  location  and  velocity 
of  the  target.  In  order  to  extend  the  relationship  between  the  output  of  the  matched  filter 
and  the  ambiguity  function  to  more  complicated  scenarios,  it  is  convenient  to  consider  the 
scatterer  distribution. 

The  scatterer  distribution,  S{T,(p),  is  the  portion  of  the  transmitted  signal  energy 
received  at  the  radar  delayed  in  time  by  an  amount  r  and  Doppler  shifted  an  amount 
<f>.  Thus,  a  target  with  location  and  velocity  which  corresponds  to  a  time  delay  tj,  and 
Doppler  shift  Wd  of  the  transmitted  signal,  causes  a  peak  in  the  scatterer  distribution  at 
{td,Wd).  The  magnitude  of  this  peak  corresponds  to  the  amount  of  energy  scattered  from 
the  target  which,  in  turn,  relates  to  the  radar  cross-section  (RCS)  of  the  target.  Figure  2.1 
shows  an  example  of  a  scatterer  distribution. 

Ideally,  the  output  of  the  matched  filter  would  be  the  scatterer  distribution,  for  it 
describes  all  the  radar  could  determine  about  the  environment.  The  transmitted  waveform 
provides  the  method  by  which  the  radar  is  able  to  measure  this  environment  and  thus  there 
exists  a  function  of  the  waveform  which  acts  as  a  filter  between  the  scatterer  distribution 
and  the  radar  receiver.  This  function  is  the  ambiguity  function  and  the  three  are  related 
by  [3]: 

M{T,4>)-=x{TA)®S{T,(j))  (11) 

where  (8  represents  2D  convolution. 

In  this  sense,  the  ambiguity  function  is  analogous  to  the  lens  of  a  camera,  which  filters 
the  true  pictme  (the  scatterer  distribution)  on  its  way  to  the  film  (the  output  of  the 
matched  filter  receiver).  By  varying  the  transmitted  waveform  it  is  possible  to  control  the 
type  of  lens  used  in  the  camera,  selecting  the  appropriate  lens  for  the  nature  of  the  objects 
being  photographed. 

It  has  been  discovered  that  there  exists  fundamental  properties  which  any  realisable 
ambiguity  function  must  obey  [2, 18].  These  are  due  to  factors  such  as  the  transmitted 
signal  possessing  positive  finite  energy  and  that  a  signal  cannot  be  defined  exactly  in  both 
time  and  frequency.  The  most  important  of  these  properties  are: 
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Example  region  of  interest 


frequenc/  shift.  4 

Corresponding  scatterer  distribution 

Figure  2.1:  Example  of  a  scatterer  distribution.  The  example  targets  have  the  following 
properties:  (a)  stationary,  medium  RCS  and  furtherest  in  range,  (b)  slowly  receding,  large 
RCS  and  closest  in  range  and  (c)  quickly  approaching,  small  RCS  and  at  medium  range. 
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•  A  Tnayimiim  of  the  magnitude  of  the  ambiguity  function  always  occurs  at  the  origin. 
That  is, 

M0,0)|>|x:(r..^)|  (12) 

•  The  magnitude  of  the  ambiguity  function  is  symmetric  about  the  line  t  =  (f).  That 

\x{-'r,-4>)\  =  \xir,(l>)\  (13) 

•  lx(0)0)l^  is  equal  to  the  amount  of  energy  in  the  transmitted  signal.  That  is, 

|x(0,0)|^=/  \x{t)\‘^dt  (14) 

Jo 

•  The  volume  under  the  siuface  of  |x(r,  ^)p  is  equal  to  the  square  of  the  energy  in  the 
transmitted  signal.  That  is, 

/  /  lx(T,0)|^dr#  = 

Jo  Jo 


These  properties  limit  the  choices  of  possible  ambiguity  functions  and  hence  possible 
radar  waveforms,  which  in  turn  affect  how  well  the  matched  filter  receiver  is  able  to 
determine  the  true  scatterer  distribution.  The  task  of  compromising  between  the  effect 
implied  by  the  ambiguity  function  and  the  necessity  to  maintain  the  above  properties  is 
the  essence  of  radar  waveform  design  and  the  subject  of  the  following  section. 


2.3  Resolution  issues 


It  can  be  seen  from  equations  (8)  and  (11)  that  in  order  to  optimally  detect  a  target 
which  produces  time  delay  and  Doppler  shift  parameters  tj,  and  respectively,  it  is 
necessary  to  maximise  the  value  of  |x(0, 0)p.  Also,  to  minimise  the  possibility  of  falsely 
detecting  targets  with  other  parameters,  it  is  necessary  to  minimise  the  value  of  |x(t,  <f>)^ 
for  (r,  ^)  7^  (0,0).  Thus  the  ideal  radar  waveform  is  one  which  produces  an  ambiguity 
function  so  that: 


Wideali-r,  <t>)? 


1  for  (t,^)  =  (0,0) 
0  otherwise 


(16) 


Unfortunately,  it  is  impossible  for  a  signal  x{t)  to  produce  this  ambiguity  function.  The 
third  and  fourth  properties  of  the  previous  section  imply  that  there  must  exist  a  constant 
volume  under  \x{'^i4J)^  which  cannot  be  confined  to  the  origin,  rendering  equation  (16) 
impossible  to  realise. 

Thus  the  ambiguity  function  must  contain  volume,  or  non-zero  magnitudes,  at  values 
of  r  and  ^  away  from  the  origin.  When  such  a  function  is  convolved  with  the  scatterer 
distribution,  three  undesirable  effects  can  occur: 


•  If  a  large  volume  of  the  ambiguity  function  is  located  near  the  origin  such  as  to 
produce  a  wide  peak  there,  it  may  prove  impossible  to  resolve  two  targets  located 
close  to  each  other  (each  target  will  produce  a  wide  peak  in  the  output  of  the  matched 
filter). 
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•  If  the  return  from  one  target  is  significantly  larger  than  the  return  from  another  and 
the  location  of  the  second  target  is  such  that  the  peak  of  the  ambiguity  function 
in  the  matched  filter  output  due  to  this  target  lies  over  a  non-zero  region  of  the 
ambiguity  function  due  to  the  first  target,  it  is  possible  that  these  non-zero  values 
are  so  large  that  they  mask  any  return  from  the  second  target. 

•  If  values  of  the  ambiguity  function  away  from  the  origin  are  large  enough,  spikes  will 
occur  in  the  output  of  the  matched  filter  at  values  of  r  and  (p  for  which  there  does 
not  correspond  a  target.  Thus  a  single  target  could  conceivably  produce  any  number 
of  false  detections. 

By  shaping  the  ambiguity  function  it  is  possible  to  reduce  the  possibiUty  of  one  or  two 
of  these  effects  from  occurring,  but  not  all  three.  If  the  central  spike  is  narrow  enough  so 
that  the  first  effect  is  avoided,  the  constant  volume  requirement  implies  that  the  second 
or  third  effect  is  likely  to  occm.  Similarly,  the  second  and  third  effects  can  be  avoided  by 
concentrating  the  volume  of  the  ambiguity  function  to  the  origin,  which  provides  the  first 
effect. 

Some  knowledge  of  the  scatterer  distribution  is  required  to  make  a  sound  judgement  on 
which  radar  waveform  to  use.  For  example,  if  it  is  known  that  all  scatterers  in  the  domain 
of  the  scatterer  distribution  have  similar  RCS,  a  reasonable  radar  waveform  would  be  one 
which  produces  an  ambiguity  function  with  a  very  narrow  spike  at  the  origin  and  the 
majority  of  the  volume  distributed  elsewhere  in  the  (r,  <j))  plane.  Similarly,  if  the  domain 
of  the  scatterer  distribution  is  confined  to  be  very  small,  a  possible  radar  waveform  is  one 
which  produces  an  ambiguity  function  which  is  zero  for  a  region  around  the  origin.  In 
this  case,  any  false  detections  due  to  the  non-zero  values  of  the  ambiguity  function  will  be 
associated  with  targets  outside  the  region  of  interest. 


2.4  Classes  of  waveforms 

Rihaczek  [4]  argues  that  any  possible  radar  waveform  could  be  considered  a  member 
of  one  of  four  classes.  Each  class  has  distinct  properties  and  ambiguity  function  shape. 

The  first  class  of  waveforms  are  those  for  which  the  product  of  their  duration,  T,  and 
bandwidth,  B,  is  approximately  umty.  The  most  basic  radar  waveform,  a  pulse  of  single 
frequency,  falls  into  this  class,  as  does  the  very  narrow,  but  wideband,  impulse  function. 
Such  waveforms  are  designated  as  Class  A  waveforms. 

The  ambiguity  function  of  a  Class  A  waveform  consists  of  a  ridge  along  either  the  line 
T  =  0  or  ^  =  0.  Thus  these  waveforms  are  able  to  provide  either  high  resolution  in  time 
delay  or  high  resolution  in  Doppler  shift,  but  not  both. 

In  contrast  to  Class  A  waveforms  are  those  signals  whose  time-bandwidth  product 
is  much  larger  than  unity.  Such  waveforms  are  often  called  pulse-compressed  signals  as 
they  are  often  generated  by  passing  a  single-frequency  pulse  through  a  dispersive  filter. 
Waveforms  of  this  type  are  divided  into  three  classes:  Class  Bl,  Class  B2  and  Class  C. 

Class  Bl  waveforms  are  irregular  or  noise- like.  Typically  they  are  aperiodic  and  uncor¬ 
related  with  themselves  over  any  more  than  a  small  delay  in  time  or  small  shift  in  Doppler. 


12 


DSTO-TR-0809 


Thus  the  ambiguity  function  of  these  waveforms  consists  of  a  spike  at  the  origin  with  the 
bulk  of  the  volume  distributed  over  a  low-level  pedestal  surrounding  the  spike.  Due  to 
their  shape,  Rihaczek  refers  to  these  ambiguity  functions  as  “thumbtack”  functions. 

Provided  the  ambiguity  function  is  normalised  so  that  |x(0, 0)|  =  1  and  that  the  volume 
contained  in  the  central  spike  is  negligible,  the  volume  of  the  pedestal  is  unity.  As  with  any 
ambiguity  function,  those  of  Class  B1  waveforms  are  confined  to  the  intervals  |t|  <  T/2 
and  \4>\  <  S/2,  so  the  average  height  of  the  pedestal  must  be  1/{TB).  This  result  is 
summarised  in  the  following  equation: 


E{|xciassBl(r, 


1  for  (r,  </>)  =  (0, 0) 
5^  otherwise 


(17) 


where  E{x}  refers  to  the  expected  value  of  x. 

It  will  be  shown  that  the  waveforms  from  transmitters  of  opportunity  are  Class  B1 
waveforms.  For  this  reason  the  ambiguity  function  of  these  signals  are  discussed  in  more 
detail  in  the  following  chapters,  including  a  more  detailed  derivation  of  the  above  equation 
in  section  4.2. 


Due  to  the  shape  of  the  ambiguity  function.  Class  B1  waveforms  typically  provide  good 
ability  to  resolve  targets  in  both  time  delay  and  Doppler  shift,  however  the  location  of 
targets  with  low  RCS  are  often  masked  by  the  pedestal  from  a  target  with  large  RCS. 

The  third  class  of  waveforms  are  Class  B2  waveforms.  Like  Class  B1  waveforms,  these 
are  aperiodic  signals  with  a  time-bandwidth  product  much  larger  than  unity,  however 
Class  B2  waveforms  are  deterministic.  They  are  typically  generated  by  modulating  a 
single  frequency  pulse  with  a  simple  function.  A  linearly  swept  FM  pulse  is  an  example 
of  a  Class  B2  waveform. 

The  ambiguity  function  of  a  Class  B2  waveform  typically  consists  of  a  ridge  orientated 
at  an  angle  to  the  r  and  (f)  axis.  Thus,  it  is  possible  with  this  waveform  to  resolve  well  all 
targets  but  those  with  a  Doppler-delay  product  which  matches  the  angle  the  ridge  is  on. 
For  these  targets,  the  ridges  produced  by  the  ambiguity  functions  will  align  and  the  exact 
location  of  the  target  on  the  (r,  (j))  would  be  difficult  to  determine. 

Class  C  waveforms  are  the  final  class  of  waveforms  and  include  all  waveforms  which  axe 
periodic.  A  periodic  signal  of  period  T  produces  an  ambiguity  function  which  is  non-zero 
only  for  values  of  (f>  which  are  multiples  of  1/T  in  the  same  way  as  the  Fourier  transform 
of  a  signal  of  period  T  is  a  set  of  lines  spaced  1/T  apart.  Also,  as  the  waveform  is  periodic, 
the  ambiguity  function  is  repeated  along  the  r  axis  over  the  same  period  T.  The  result  is 
a  set  of  spikes  spaced  1/T  apart  along  the  0  axis  and  spaced  at  most  T  along  the  r  axis. 
This  shape  is  often  referred  to  as  a  “bed  of  nails”. 

Class  C  waveforms  thus  provide  good  resolution  in  both  time  delay  and  Doppler  shift 
and  are  able  to  resolve  targets  of  greatly  varying  RCS,  however  the  matched  filter  output 
can  contain  multiple  spikes  for  any  single  target,  thus  producing  an  ambiguity. 

Example  ambiguity  functions  of  the  four  classes  of  waveforms  are  shown  in  figure  2.2. 
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Class  A  Class  B1 


Figure  2.2:  Typical  ambiguity  functions  of  the  four  classes  of  waveforms. 


2.5  Implementation  issues 

The  previous  sections  provide  detail  on  how  the  shape  of  the  waveform  affects  how 
well  the  radar  is  able  to  resolve  targets.  In  practice  however,  much  of  the  processing  is 
implemented  on  digital  computers  and  this  produces  further  limits  in  target  resolution. 

Equation  (6)  describes  how  the  matched  filter  is  implemented  on  a  computer.  The 
signals  Xr[n]  and  Xo[n]  are  the  received  and  original  radar  waveforms  respectively,  sampled 
at  the  discrete  points  in  time  nTg,  where  Tg  is  the  sampling  period.  Prom  this  equation  it 
can  be  seen  that  it  does  not  make  sense  to  evaluate  M (t,  for  non-integer  values  of  t,  as 
Xo[n]  is  not  defined  for  non-integer  values  of  n.  This  corresponds  to  being  able  to  compare 
the  received  signal  with  the  transmitted  one  only  delayed  in  time  by  multiples  of  T^. 

Also,  the  fact  that  the  received  signal  has  finite  dmation,  T  =  NTs  s,  implies  a  degree  of 
smoothness  in  its  spectrum.  Prom  Fourier  theory  is  it  known  that  such  a  spectrum  contains 
independent  information  only  at  intervals  of  1/T,  so  it  makes  sense  to  evaluate  M{T,(f>) 
only  for  values  of  (f)  which  correspond  to  shifts  in  frequency  of  multiplies  of  l/{NTs)  Hz. 

Of  course,  it  is  possible  to  evaluate  M(r,  if)  for  values  of  r  ^  kTs  and  for  values  of 
(f  ^  k/{NTs)  (where  k  is  an  integer),  however  the  result  is  merely  an  interpolation  of  the 
values  of  M (r,  (f)  which  he  on  this  evenly  spaced  grid. 

This  implies  that  it  is  not  possible  to  resolve  in  range  targets  which  produce  a  signal 
return  delayed  less  than  Tg  s  apart,  nor  is  it  possible  to  resolve  in  velocity  targets  which 
produce  a  signal  return  shifted  in  Doppler  less  than  l/^NTg)  Hz.  This  corresponds  to 
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targets  no  less  than  cTj  m  apart  and  travelling  with  velocities  different  by  no  less  than 
c/{NTsfc),  where  c  is  the  speed  of  light  and  fc  is  the  carrier  frequency  of  the  waveform. 

To  achieve  better  range  resolution  it  is  required  to  decrease  Tg,  the  sample  period,  or 
rather,  increase  the  frequency  at  which  the  radar  waveform  is  sampled.  Of  course,  the  radar 
waveform  must  have  energy  at  these  higher  frequencies,  which  corresponds  to  Woodward’s 
claim  that  range  resolution  is  a  function  of  the  bandwidth  of  the  radar  waveform.  To 
achieve  better  resolution  in  velocity  it  is  required  to  increase  T  =  NTg,  the  length  of  radar 
waveform  used  by  the  processor.  By  increasing  the  sampling  frequency  and  length  of  Xr\p\ 
and  Xo[n],  the  grid  upon  which  M(r,  4>)  can  be  evaluated  becomes  finer. 

Increasing  the  sampling  frequency  or  integration  time  assists  the  ability  to  resolve 
targets  in  a  second  way.  Usually  the  received  signal,  a;r[n],  is  corrupted  with  noise  and 
M(r,  contains  a  random  component.  Increasing  the  sampling  frequency  or  integration 
time  implies  more  samples  of  SrW  are  being  used  and  the  variance  of  M(t,  decreases. 
If  interpolation  is  used  to  gain  a  more  accurate  estimate  of  the  time  delay  or  Doppler  shift 
of  a  target,  this  estimate  is  improved. 

It  should  be  mentioned  again  that  target  resolution  is  a  function  of  both  the  shape  of 
the  waveform  and  the  granularity  of  this  digital  grid.  It  is  futile  increasing  the  sample 
rate  or  integration  time  to  produce  a  grid  so  fine  that  the  width  of  the  peak  at  the  origin 
of  the  ambiguity  function  encompasses  hundreds  of  (r,  (f>)  cells. 
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Chapter  3 


Non-cooperative  Passive  Radar 


This  chapter  describes  the  basic  elements  of  a  radar  which  uses  the  transmissions  already 
present  in  the  region  of  interest  as  the  radar  waveform.  The  transmitters  being  utilised  are 
often  referred  to  as  non-cooperative  since  their  existence  and  characteristics  are  dictated 
by  their  own  missions  rather  than  by  those  of  the  radar. 

Section  3.1  outlines  the  proposed  radar  configuration  which  utilises  these  transmissions. 
It  defines  the  system  which,  throughout  the  remainder  of  the  report,  will  be  investigated 
for  feasibility. 

This  configuration  is  clearly  bistatic,  as  the  transmitter  is  located  away  from  the  re¬ 
ceiver.  Section  3.2  discusses  some  essential  principles  of  bistatic  radar,  specifically  how 
the  parameters  gained  from  the  output  of  the  matched  filter  can  be  used  to  determine  the 
actual  location  and  velocity  of  the  target. 

Section  3.3  describes  what  transmissions  of  opportunity  exist.  Some  initial  character¬ 
istics  of  these  waveforms  are  also  given,  however  a  more  detailed  analysis  is  left  for  the 
following  chapter. 


3.1  Radar  configurations 

The  theory  of  the  previous  chapter  explains  that  two  signals  are  required  to  perform 
matched  filter  reception:  a  copy  of  the  transmitted  signal  and  the  received  signal.  An 
obvious  potential  problem  of  non-cooperative  passive  radar  is  that  there  is  no  copy  of  the 
transmitted  signal  at  the  receiver. 

To  overcome  this,  it  first  must  be  noted  that  the  transmission  of  opportunity  typically 
arrives  at  the  radar  receiver  from  a  specific  direction.  Usually,  this  is  the  direction  from 
the  receiver  directly  to  the  transmitter,  although  under  certain  propagation  environments 
this  may  not  be  so. 


3.1.1  Two  receiver  configuration 

The  first  proposed  radar  system  uses  two  receivers.  One  receiver  (and  its  corresponding 
antenna)  is  fixed  to  collect  the  signal  arriving  at  the  radar  travelling  directly  from  the 
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Transmitter 


transmitter  of  opportunity.  This  signal  will  be  used  as  the  copy  of  the  transmitted  signal 
and  will  be  referred  to  as  Xo{t).  The  second  receiver  and  antenna  is  configinred  to  gather 
energy  from  the  bearing  along  which  the  detection  of  targets  is  being  attempted.  The 
signal  received  from  this  direction  will  be  referred  to  as  Xr{i). 

It  is  important  that  Xo[t)  contain  as  much  of  the  direct  signal  and  as  little  of  any  other 
signal  as  possible,  and  that  Xr  (t)  contain  as  much  of  the  signal  arriving  from  the  direction  of 
interest  and  as  little  from  any  other  direction  as  possible.  The  responsibility  for  performing 
this  task  falls  to  the  directivity  of  the  antennas  and  the  ability  to  attenuate  signals  arriving 
from  directions  other  than  the  direction  of  interest.  This  concept  is  illustrated  in  figure  3.1. 

Where  this  presents  the  largest  problem  is  at  the  receiver  scanning  for  target  scattering. 
Typically  the  energy  scattered  from  a  target  is  less  than  that  in  the  direct  signal,  so  unless 
the  antenna  linked  to  this  receiver  attenuates  sufficiently  in  the  direction  towards  the 
transmitter  of  opportunity,  Xr{t)  could  conceivably  contain  more  energy  from  the  direct 
signal  than  from  the  scattered  signal. 

One  method  for  controlling  this  energy  leakage  is  to  use  an  array  of  antennas  and 


18 


DSTO-TR-0809 


Direct  signal  / 
from  transmitter  / 


Beams  formed 
in  the  direction  of 
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Array  of  receivers 


Figure  3.2:  Electronically  beam-steering  with  null  placement  to  both  scan  in  multiple  di¬ 
rections  simultaneously  and  minimise  the  energy  leaking  from  the  direct  signal. 


electronically  form  a  beam  in  the  the  direction  of  interest  so  that  a  null  is  placed  in  the 
bearing  along  which  the  direct  signal  is  arriving.  This  method  also  provides  the  ability  to 
form  multiple  beams  and  effectively  look  in  a  number  of  different  directions  simultaneously. 
Figure  3.2  illustrates  this  technique.  The  study  of  cirray  processing  and  beam-forming  is 
a  matiue  field  and  many  procedures  have  been  developed  for  performing  this  type  of 
processing  [19,20]. 


3.1.2  Single  receiver  configuration 

The  second  passive  radar  configuration  under  investigation  uses  only  a  single  receiver. 
In  this  configuration  the  received  signal  would  contain  energy  both  firom  the  direct  signal, 
Xo{t),  and  from  the  scattered  signal,  Xr{t).  It  will  be  shown  that  for  this  configuration 
to  be  effective,  the  gain  of  the  antenna  in  the  direction  of  interest  needs  to  be  far  greater 
than  that  in  the  direction  towards  the  transmitter  of  opportunity.  The  single  receiver 
configiuation  is  illustrated  in  figure  3.3. 

In  this  configuration,  the  single  received  signal  is  used  both  as  Xo{t)  and  Xr{t)  when 
performing  matched  filtering. 


3.1.3  Other  configurations 

For  the  radar  configurations  of  the  previous  paragraphs  to  be  able  to  resolve  the  lo¬ 
cation  or  velocity  of  the  target  it  is  necessary  to  accurately  measure  the  bearing  of  the 
target,  6.  To  accomplish  this  for  the  values  of  carrier  frequency  proposed^,  the  necessary 
antennas  usually  need  to  be  either  very  large  or  very  complex. 

^Details  of  the  proposed  transmissions  of  opportunity,  including  their  carrier  frequencies,  can  be  found 
in  section  3.3 
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Transmitter 


Figure  3.3:  Configuration  of  the  passive  radar  using  a  single  receiver. 
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Transmitter 


Figure  3. 4:  Configuration  of  the  passive  radar  using  a  two  receivers  located  apart. 

A  third  proposed  passive  radar  configuration  involves  two  receivers  located  far  apart, 
each  using  antennas  which  provide  poor  resolution  in  direction,  but  are  still  capable  of 
heavily  attenuating  the  signal  arriving  from  the  direction  of  the  transmitter  of  opportunity. 
Such  antennas  are  typically  small  and  inexpensive  to  produce.  The  proposed  configuration 
is  illustrated  in  figure  3.4. 

In  this  configuration,  each  receiver  performs  its  own  matched  filtering  using  the  single 
signal  in  the  exactly  the  same  manner  as  the  receiver  in  the  single  receiver  configuration. 
In  this  system  however,  direction  of  the  scattered  signal,  $,  is  unknown  and  the  location 
and  velocity  of  the  target  is  resolved  using  the  second  receiver. 

Another  possible  configuration  involves  utilising  more  than  one  transmitter  of  oppor¬ 
tunity.  This  multi-static  arrangement  would  consist  of  a  number  of  receivers  and  corre¬ 
sponding  antenna,  each  configured  to  one  of  the  three  previously  proposed  configurations. 
Each  receiver  would  perform  its  own  matched  filtering,  the  results  of  which  would  then  be 
combined  in  some  manner  to  provide  a  single  estimate  of  the  number  of  targets  present, 
their  location  and  their  velocity. 
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3.2  Bistatic  considerations 


Chapter  2  showed  that  matched  filter  processing  provides  information  on  how  much 
the  signal  scattered  from  a  target  is  delayed  in  time  and  shifted  in  Doppler.  The  conversion 
between  the  signal  delay,  ta,  and  frequency  shift,  Wd,  produced  by  a  target  and  the  target’s 
actual  location  and  velocity  depends  on  the  configuration  of  the  radar.  This  section 
examines  this  conversion  for  the  systems  proposed  in  the  previous  section. 


Figures  3.1  and  3.3  illustrate  the  dual  and  single  receiver  configurations  proposed  to 
utilise  transmissions  of  opportunity.  It  is  assumed  that  9  and  tq  are  known,  so  in  order  to 
determine  the  location  of  the  target,  only  the  the  distance  r2  is  required,  which  is  given 
by: 


ctdjctd  +  2ro) 

2(ctrf +  ro(l  -COS0)) 


(18) 


where  c  is  the  speed  of  light. 


The  velocity  of  the  target,  v,  and  direction  the  target  is  travelling,  are  given  by: 


and 


where 


V  =  Vri+  v‘^2  +  2Url^^r2  COs(0') 


=  arcsin 


^sin(0') 

V 


(19) 

(20) 


Vr\  =  target’s  velocity  in  the  direction  of  ri 

r2  f  cwd  .  dd  .  ^  ^ 

=  - ; - hr-y +ror2— sin0  -  rocos0 

n+T’2V27r/c  dt  J 

Vr2  =  target’s  velocity  in  the  direction  of  r2 

-ri  f  cwd  ,  dd  .  „ 

-  — ; — hr^+’"o’'237sm0-rocos0 

ri  +  r-2V27r/c  dt  J 

o'  =  angle  between  ri  and  r2 
ro 


=  arcsin 


ri 


sin0 

Ln 

ro  -  r2  +  ctd 


(21) 

(22) 

(23) 

(24) 


and  fc  is  the  carrier  frequency  of  the  radar  transmission.  The  angle  9^  is  measured  relative 
to  the  direction  towards  the  target.  For  example,  if  the  target  is  travelling  directly  away 
from  the  radar,  =  0.  The  derivation  of  equations  (18)  through  (24)  can  be  found  in 
appendix  A. 

Note  that  unlike  monostatic  radar,  the  ability  to  measure  the  range  of  the  target,  r2, 
is  a  function  of  the  bearing  to  it,  9,  which  may  prove  difficult  to  measure  accurately.  Also, 
the  ability  to  accurately  measure  the  target’s  velocity  depends  on  the  ability  to  accurately 
measure  its  range  (equations  (19)  and  (20)  are  functions  of  r2).  This  imphes  the  ability 
to  resolve  well  both  the  delay  in  time,  td,  and  the  Doppler  shift,  Wd,  caused  by  the  target. 

Often  these  inaccuracies  can  be  minimised  by  using  a  number  of  observations  and  an 
adaptive  or  block-based  processing  technique  with  the  objective  of  optimally  estimating 


22 


DSTO-TR-0809 


the  target  location  and  velocity  against  a  defined  criteria  (such  as  least  square  error).  The 
radar  proposed  by  Howland  [13]  used  such  a  method  for  resolving  target  location  using 
a  number  of  observations  of  Doppler  shift,  Wd,  as  the  waveform  used  was  such  that  time 
delay,  td,  was  not  able  to  be  measured  accurately. 

Techniques  for  optimal  parameter  estimation  are  well  developed  and  the  reader  is 
referred  to  Anderson  and  Moore  [21]  for  more  information. 

Up  to  this  point,  it  has  been  assumed  that  the  location  of  the  transmitter  of  opportunity 
is  known.  This  is  usually  the  case,  as  the  location  of  radio  and  television  broadcast  towers 
are  readily  attainable.  If  however,  these  details  are  not  available,  they  are  simply  calculated 
using  two  receivers  located  apart.  Each  receiver  scans  for  the  direction  from  which  the 
strongest  signal  arrives.  Typically,  this  is  the  direct  signal  and  lines  drawn  fi-om  each 
receiver  outwards  in  this  direction  will  bisect  at  the  location  of  the  transmitter. 

For  a  more  detailed  discussion  on  the  implications  of  bistatic  radar,  including  equations 
for  resolving  target  location  and  velocity  for  more  complex  scenarios  than  those  described 
above,  the  reader  is  referred  to  Skolnik  [22]  or  Willis  [23]. 

3.3  Waveforms  of  opportunity 

The  electromagetic  spectrum  is  abundant  in  signals  which  arise  fi'om  numerous  sources 
such  as  commercial  broadcasting  stations,  private  and  government  radar  transmitters  and 
various  communications  systems.  Of  these,  this  report  proposes  to  investigate  the  following 
transmissions  for  their  suitability  to  passive  radar^: 

•  AM  radio 

—  Commercial  broadcast  radio: 

Signals  of  up  to  20  kHz  bandwidth  occurring  at  firequencies  between  530  kHz 
and  1605  kHz  (wavelengths  of  187  m  to  566  m) 

—  Medium  wave  and  short  wave  radio: 

Signals  of  up  to  15  kHz  bandwidth  occmring  at  various  frequencies  between 
1605  kHz  and  26.9  MHz  (wavelengths  of  11  m  to  187  m) 

•  FM  radio 

—  Commercial  broadcast  radio: 

Signals  of  up  to  75  kHz  bandwidth  occurring  at  frequencies  between  88  MHz 
and  108  MHz  (wavelengths  of  2.8  m  to  3.4  m) 

—  Other  frequency  modulated  communications  originating  firom  a  variety  of 
somces  such  as  land,  sea  and  air  mobile  communications  equipment  and  satellite 
transmissions 

•  Television  broadcasts 

Signals  of  up  to  6  MHz  bandwidth  occurring  at  the  following  frequencies: 

^Details  of  frequencies  and  bandwidths  used  in  this  list  come  from  [24]. 
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-  Channels  0-4: 

46.25  MHz  to  95.25  MHz  (wavelengths  3.1  m  to  6.5  m) 

—  Channels  6-11: 

175.25  MHz  to  216.25  MHz  (wavelengths  1.4  m  to  1.7  m) 

-  UHF  channels: 

527.25  MHz  to  814.25  MHz  (wavelengths  0.4  m  to  0.6  m) 

The  characteristic  which  distinguishes  these  transmissions  is  that  they  are  all  random 
aperiodic  signals.  They  are  generated  from  random  information  such  as  speech,  music 
and  video,  and  modulated  using  one  of  three  schemes:  frequency  modulation,  amplitude 
modulation  or  television’s  vestigial  sideband  modulation. 

This  random  information  is  typically  uncorrelated  with  itself  over  any  more  than  a 
small  delay  in  time  or  shift  in  frequency.  Plain  speech,  for  example,  is  uncorrelated  with 
a  copy  of  itself  delayed  any  more  than  10  ms  [25]. 

The  modulation  applied  to  these  signals  modifies  this  structure.  For  the  three  schemes 
mentioned  above,  the  modulation  causes  the  energy  in  the  original  signal  to  be  distributed 
over  a  wider  range  of  frequencies  and  in  the  case  of  frequency  modulation,  to  be  distributed 
more  evenly  across  these  frequencies.  Thus,  provided  the  pass-band  of  the  receiver  is 
smaller  than  the  bandwidth  of  the  signal  produced  by  this  modulation,  the  spectrum  of 
the  received  signal  contains  energy  distributed  more  evenly  across  all  retained  frequencies, 
and  thus  is  uncorrelated  with  itself  delayed  over  even  smaller  shifts  in  time  and  frequency, 
than  the  original  signal. 

This  concludes  that  the  received  signal  of  opportunity  is  not  unlike  bandlimited  white 
noise.  This  assumption  is  strengthened  when  the  receiver  scans  for  long  periods  of  time 
and  any  underlying  structure  of  the  signal  due  to  the  nature  of  the  transmission  is  averaged 
out. 

This  characteristic  is  most  apparent  in  the  ambiguity  function  of  these  signals.  Sec¬ 
tion  2.4  showed  that  the  ambiguity  function  of  a  random  signal  contains  a  spike  at 
(r,  =  (0,0)  surrounded  by  a  constant  pedestal  of  energy.  Figure  3.5  shows  the  am¬ 
biguity  function  of  an  example  transmission  of  opportimity.  The  signal  used  to  generate 
this  figure  was  taken  from  taken  from  a  local  FM  radio  station. 

The  similarity  between  the  ambiguity  function  for  Class  B1  waveforms,  shown  in  fig¬ 
ure  2.2,  and  figure  3.5  confirm  that  waveforms  from  the  transmitters  of  opportunity  being 
investigated  fall  into  Rihaczek’s  second  class  of  waveforms. 

The  similarity  between  the  transmission  of  opportunity  and  white  noise  breaks  down  as 
the  bandwidth  of  the  receiver  approaches  the  bandwidth  of  the  transmission  of  opportunity. 
The  result  is  that  the  amount  of  time  and  frequency  shift  required  before  the  signal  is 
uncorrelated  with  itself  increases  and  the  spike  in  the  ambiguity  function  at  (r,  (f))  =  (0, 0) 
becomes  wider. 

As  the  ability  to  resolve  targets  is  a  function  of  the  width  of  the  spike  at  the  origin  of 
the  ambiguity  function,  it  is  desired  to  maintain  this  similarity  to  white  noise,  however  it 
is  also  desired  to  use  a  signal  of  large  bandwidth  so  that  the  width  of  a  cell  along  the  (p 
axis  is  kept  to  a  minimum  (as  explained  in  section  2.5). 
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Figure  3.5:  Ambiguity  function  from  an  example  transmission  of  opportunity.  The  3  s 
long  example  signal  was  captured  from  a  local  FM  radio  broadcast  using  a  receiver  which 
had  4  kHz  bandwidth. 

Further  discussion  on  the  characteristics  of  the  transmissions  of  opportunity  requires 
knowledge  of  the  individual  modulation  schemes.  The  remainder  of  this  section  discusses 
the  effect  of  these  three  schemes. 


3.3.1  AM  radio 

Amplitude  modulation  produces  the  signal  with  the  smallest  bandwidth  of  the  three 
types  of  transmissions  being  investigated.  It  also  produces  the  signal  which  contains  the 
most  spectral  variance  over  the  same  amount  of  bandwidth.  Thus,  it  is  least  capable  of 
being  modelled  by  white  noise  and  would  most  likely  make  for  the  worst  radar  waveform. 

More  than  half  of  the  energy  of  an  AM  transmission  is  contained  in  the  carrier.  Most  of 
the  remaining  energy  resides  between  the  frequencies  /c+100  and  /c+3000  Hz  and  between 
fc  —  100  and  fc  —  3000  Hz,  where  /c  is  the  frequency  of  the  carrier.  This  distribution  of 
energy  results  in  the  received  signal  being  not  as  white  as  those  produced  by  the  other 
modulation  schemes  and  the  corresponding  ambiguity  function  not  containing  as  distinct 
a  peak  at  the  origin. 

An  example  of  the  ambiguity  function  of  a  typical  AM  radio  transmission  is  shown  in 
figure  3.6.  In  this  example,  the  modulation  has  not  removed  the  underlying  harmonics 
and  structure  contained  in  the  original  signal  (which  was  plain  speech)  and  the  relatively 
large  bandwidth  has  been  able  to  capture  this.  This  effect  is  evident  in  the  large  peaks 
which  are  located  away  from  the  origin  in  the  ambiguity  function. 

Obviously,  this  effect  is  undesired  in  a  radar  waveform.  As  mentioned  in  section  2.3, 
these  peaks  could  mask  smaller  targets  or  produce  false  detections. 
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(b) 


Figure  3.6:  Image  (a)  and  mesh  (b)  view  of  the  ambiguity  function  of  an  example  AM 
radio  transmission.  The  example  signal  was  0.7  s  long  and  supported  a  4  kHz  bandwidth. 
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Figure  3.7:  Image  (a)  and  mesh  (b)  view  of  the  ambiguity  function  of  an  example  AM  radio 
transmission.  The  example  signal  was  0.7  s  long  and  was  limited  to  a  1  kHz  bandwidth. 


One  method  to  improve  the  shape  of  this  ambiguity  function  is  to  take  a  longer  time 
sample,  which  would  result  in  the  harmonics  or  structure  to  be  averaged  out. 

Another  method  of  attempting  to  shape  the  waveform  into  a  more  usable  form  is  to  use 
a  smaller  bandwidth.  This  has  the  effect  of  moving  the  energy  of  the  ambiguity  function 
into  the  centre  however  also  causing  this  centre  peak  to  stretch  along  the  r  axis.  An 
example  of  the  ambiguity  function  of  the  same  signal  used  in  figure  3.6,  passed  through  a 
filter  to  limit  its  bandwidth  to  1  kHz,  is  shown  in  figure  3.7. 

The  extreme  of  this  technique  is  a  receiver  whose  bandwidth  is  so  narrow  that  Xo{t) 
effectively  consists  of  a  single  frequency  and  the  ambiguity  function  looks  like  a  ridge  along 
the  T  axis,  like  class  A  waveforms  from  section  2.4. 
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3.3.2  FM  radio 

Frequency  modulation  produces  a  signal  much  closer  to  the  desired  bandlimited  white 
noise  model  than  does  amplitude  modulation,  and  over  a  wider  bandwidth.  The  spectrum 
of  an  FM  radio  transmission  extends  for  75  kHz  with  the  energy  being  well  distributed 
across  most  of  this  bandwidth. 

Figures  3.8,  3.9  and  3.10  show  the  ambiguity  functions  of  example  FM  radio  transmis¬ 
sions.  In  each  figure,  a  different  signal  was  being  modulated  and  transmitted:  figure  3.8 
shows  the  ambiguity  function  of  an  FM  transmission  of  rock  music,  figme  3.9  of  classical 
music  and  figure  3.10  of  plain  speech.  It  can  be  seen  that  there  is  not  a  lot  of  variation  be¬ 
tween  the  three  ambiguity  functions,  although  fi:equency  modulated  plain  speech  appears 
to  provide  the  finest  resolution  of  Doppler  shift  (figure  3.10  contains  the  narrowest  ridge). 


3.3.3  Television  broadcasts 

The  modulation  scheme  used  for  television  broadcasts  is  slightly  more  complex  than 
that  used  for  AM  or  FM  radio^.  Three  signals  are  coded  together  over  the  6  MHz  band¬ 
width  used  by  the  transmission:  a  chrominance  signal,  a  luminance  signal  and  the  sound 
track.  Mixed  with  this  however,  are  a  number  of  pulses  used  by  the  receiving  equipment 
to  determine  various  synchronising  information.  The  most  dominant  of  these  is  the  line 
sync  pulse. 

The  line  sync  pulse  is  4.7  ^J,s  long,  occurs  every  64  fj,s  and  has  an  amplitude  at  least  132 
percent  of  any  other  information  in  the  television  broadcast.  For  this  reason  it  is  expected 
that  the  ambiguity  function  of  a  television  broadcast  contains  narrow  ridges  extending 
±l/4.7/xs  =  ±213  kHz  in  the  4>  dimension  repeated  every  64  //s  in  the  r  dimension. 

Thus  the  return  from  a  single  target  produces  any  number  of  returns,  each  separated 
along  the  r  axis  by  64  fxs.  If  the  region  of  interest  is  small  enough,  this  ambiguity  may 
cause  no  problem. 

This  effect  is  also  removed  if  the  bandwidth  of  the  receiver  is  too  small  to  detect  this 
15,625  times-a-second  pulse.  When  received  using  filter  bandwidths  smaller  than  15  kHz, 
the  transmission  appears  to  follow  closely  the  desired  white  noise  model.  However,  as 
mentioned  in  chapter  2,  there  are  disadvantages  of  using  smaller  bandwidths.  Figure  3.11 
shows  the  ambiguity  function  from  a  typical  television  broadcast  taken  using  a  receiver 
which  had  a  4  kHz  bandwidth. 


^Details  of  the 


modulation  scheme  used 


by  television  comes  from 


[26]. 
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Figure  S.8:  Image  (a)  and  mesh  (b)  view  of  the  ambiguity  function  of  an  example  FM 
transmission  of  rock  music.  The  example  signal  was  2.7  s  long  and  supported  a  12  kHz 
bandwidth. 
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(a) 


Doppler  shift,  i|>  (Hz) 


(b) 


Figure  3.9:  Image  (a)  and  mesh  (b)  view  of  the  ambiguity  function  of  an  example  FM 
transmission  of  classical  music.  The  example  signal  was  2.7  s  long  and  supported  a  12  kHz 
bandwidth. 
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Figure  3.10:  Image  (a)  and  mesh  (b)  view  of  the  ambiguity  function  of  an  example  FM 
transmission  of  speech.  The  example  signal  was  2.7  s  long  and  supported  a  12  kHz  band¬ 
width. 
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Figure  3.11:  Ambiguity  function  from  an  example  television  broadcast.  The  3  s  long 
example  signal  was  captured  using  a  receiver  which  had  4  kHz  bandwidth. 
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Chapter  4 

Analysis  of  Waveforms  of 
Opportunity 


This  chapter  combines  the  radar  system  and  particular  waveforms  of  opportunity  proposed 
in  the  previous  chapter  with  the  theory  given  in  chapter  2.  The  approach  taken  is  to 
calculate  statistical  properties  of  the  output  of  the  matched  filter  receiver,  specifically  the 
mean  and  variance,  and  determine  whether  it  would  be  possible  to  accurately  resolve  the 
location  and  velocity  of  targets  from  it. 

Section  4.1  states  more  formally  the  white  noise  model  of  the  waveforms  of  opportunity, 
which  is  used  throughout  the  remainder  of  the  chapter.  Section  4.2  applies  this  model  to 
derive  statistical  properties  of  the  ambiguity  function,  y(T,  ^),  of  these  waveforms.  These 
results  are  called  upon  in  sections  4.3  and  4.4  which  calculate  similar  properties  for  the 
output  of  the  matched  filter  and  square  law  device,  |M(r,(/»)p. 

The  mean  and  variance  of  |M(t,  provide  insight  into  the  effectiveness  of  the  trans¬ 
missions  of  opportunity  as  radar  waveforms.  Section  4.5  discusses  these  results  and  outlines 
how  other  parameters  of  the  radar  can  be  modified  to  improve  its  performance. 

The  mean  and  variance  of  the  output  of  the  matched  filter  receiver  also  provide  a 
method  of  calculating  a  threshold  against  which  target  detection  can  be  performed.  This 
is  the  subject  of  section  4.6,  which  completes  the  algorithm  for  converting  the  received 
signals  into  a  list  of  target  locations  and  velocities. 


4.1  Model  of  waveforms  of  opportunity 

As  mentioned  in  section  3.3,  the  waveforms  received  from  the  proposed  transmitters 
of  opportunity  are  random,  aperiodic  and  uncorrelated  with  themselves  when  delayed  any 
more  than  a  small  amount  in  time  or  shifted  any  more  than  a  small  amount  in  frequency. 
An  approximate  model  for  such  a  signal  is  bandlimited  Gaussian  white  noise. 

For  the  analysis  of  this  chapter,  the  sampled  version  of  the  waveform  broadcasted  from 
the  transmitter  of  opportunity,  a:[n],  is  modelled  as  a  Gaussian  white  noise  process  with 
zero  mean  and  variance  cr^.  That  is: 

E{a:[n]}  =  0  (25) 
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(26) 

(27) 


E{x^[n]} 

E{a;[n]a;*[n  + 


Ofor  {T,<f>)  ^  (0,0) 


4.2  Analysis  of  the  ambiguity  function 

It  is  desired  to  analyse  (/>)  when  the  signal  of  interest  is  the  transmission  of  oppor- 
tumty,  a;[n],  as  described  in  the  previous  section.  Prom  equation  (10): 

1  ^ 

E{x(r,</.)}  =  ;^5;E{a;[n]a:*[n  +  T]e-*^} 

n=l 

=  f  for  (t,(^)  =  (0,0) 

(  0  otherwise  ^  ^ 

where  N  is  the  number  of  samples  in  the  signal.  It  can  be  further  shown  that: 


E{|x(r,<^)P}  = 


(1  +  f  Vx  for  (r,<^)  =  (0,0) 
^(7^  otherwise 


(29) 


This  equation  gives  an  indication  to  the  shape  of  the  the  ambiguity  function  of  a 
transmission  of  opportunity.  In  order  to  compare  this  function  to  the  predicted  function 
of  equation  (17)  and  the  example  observed  functions  of  section  3.3,  this  equation  needs  to 
be  normalised  so  that  |x(0, 0)p  =  1.  Normalising  equation  (29)  provides  a  spike  of  unit 
height  at  |x(0,0)p  surrounded  by  a  pedestal  of  value: 


E{|x(r,^6)P} 


N^x  ^  i-  ^ 

(l  +  i)a4  l  +  N'^N 


(30) 


for  large  N  and  (t,  ^)  7^  (0, 0). 

Given  that  the  received  signal  was  sampled  at  the  Nyquist  frequency  necessary  to 
maintain  the  bandlimited  white  noise  model,  the  time  between  consecutive  samples  is 
f/B,  where  B  is  the  bandwidth  of  the  signal.  If  the  length  of  the  signal  is  T,  the  number 
of  samples  is  iV  =  TB.  Thus  the  height  of  the  expected  pedestal  around  the  central  spike 
in  the  ambiguity  function  is  1/(T5),  confirming  equation  (17). 

These  expected  values  for  the  ambiguity  function  match  well  with  the  ambiguity  func¬ 
tions  of  the  actual  transmissions  of  opportunity  presented  in  section  3.3.  Figure  3.5  shows 
a  central  spike  surrounded  by  a  pedestal  whose  average  value  lies  around  —40  dB.  Given 
that  the  signal  used  to  generate  this  figure  had  length  T  =  3  s  and  bandwidth  J5  =  4  kHz, 
the  height  of  this  pedestal  is  expected  to  be  1/(TB)  =  8.3  x  10“®  =  —40.8  dB. 

Similarly,  the  signal  used  to  generate  figure  3.11  was  also  of  length  T  =  3  s  and 
bandwidth  B  =  A  kHz,  implying  a  constant  level  pedestal  in  the  ambiguity  function  also 
of  —40.8  dB,  which  is  apparent  in  the  figure. 
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Figures  3.8,  3.9  and  3.10  were  each  generated  using  signals  of  length  T  =  2.7  s  and 
bandwidth  JB  =  12  kHz.  Thus,  it  is  expected  that  these  ambiguity  functions  contain  a 
central  spike  smrounded  by  a  pedestal  of  level  1/(T5)  =  3.1  x  10~®  =  —45.1  dB.  Again, 
this  agrees  with  what  is  seen  in  these  images. 

Up  to  this  point,  only  the  mean  value  of  the  ambiguity  function  has  been  considered. 
For  further  insight  into  the  statistical  properties  of  this  random  variable,  it  is  worth  cal¬ 
culating  its  variance.  This  also  will  prove  valuable  when  considering  a  detection  threshold 
for  the  output  of  the  matched  filter. 

It  can  be  shown  that: 


e{|x(t,^)|‘}  = 


■{ 


(1  4-  ^  -I-  jl?  -1-  for  (r,  0)  —  (0, 0) 


otherwise 


So  that: 


(31) 


VAR{|x(T,^i)|2}  = 


I 


iw  W  w)^x  for  (r,  <f))  =  (0, 0) 
otherwise 


(32) 


The  calculations  used  to  derive  these  results  can  be  found  in  Appendix  A.3. 

These  equations  provide  first  insights  as  to  how  other  parameters  of  the  radar  can  be 
modified  to  improve  the  ability  to  resolve  targets.  Prom  equation  (30)  it  can  be  seen  that 
by  increasing  the  number  of  samples  in  the  received  signal,  N,  the  level  of  the  pedestal  in 
the  ambiguity  function  decreases.  This  reduces  the  possibility  that  the  returns  firom  one 
target  mask  the  retmns  from  a  second  target  with  lower  RCS. 

These  implications  axe  discussed  further  in  section  4.6. 


4.3  Result  of  the  matched  filter  (single  target) 

Consider  the  proposed  passive  radar  configuration  of  section  3.1  which  utilises  two 
receivers.  If  there  exists  a  single  target  in  the  region  of  interest  the  two  signals  received 
are: 


Xo[n]  =  aox[n]  +  0!ia;[n  —  ti]e  (33) 

Xr[n]  =  /3ox[n]  + /3ix[n  —  (34) 

where  ti  and  wi  are  the  delay  in  time  and  shift  in  Doppler  of  the  transmitted  signal  due 
to  the  location  and  velocity  of  the  target. 

Note  that  noise  is  not  included  in  the  model.  In  the  following  sections  it  will  be 
shown  that  target  detection  is  very  dependent  on  interference  from  the  individual  signal 
components  of  the  received  signals  (clutter)  and  it  is  believed  that  the  effect  of  noise  is 
negligible  compared  with  this  problem. 
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The  coefficients  Ok  and  account  for  the  attenuation  of  the  antenna  in  the  direction 
of  arrival  of  the  two  components  of  the  received  signals.  As  mentioned  in  section  3.1, 
the  antenna  linked  to  the  receiver  collecting  Xo[n]  should  be  arranged  so  that  all  signals 
but  those  arriving  directly  from  the  transmitter  of  opportunity  are  attenuated,  and  the 
antenna  linked  to  the  receiver  collecting  Xr\n\  be  arranged  so  that  all  signals  but  those 
arriving  from  the  bearing  being  investigated  for  targets  are  attenuated. 

The  coefficients  cxk  and  Ph  also  account  for  the  loss  and  phase  shift  which  occurs  when 
the  signal  scatters  from  the  target  and  the  loss  due  to  the  length  of  the  path  travelled  by 
the  signal  from  the  transmitter  to  the  receiver. 

Thus  it  is  expected  that  ao  is  large,  j3i  is  smaller  and  ai  and  /Sq  are  close  to  zero. 
Note  that  due  to  their  dependence  on  the  location  of  the  targets,  ak  and  are  random 
variables. 

The  result  of  the  matched  filter  receiver  is  found  by  applying  the  above  signals  to 
equation  (5): 


1  ^ 

M{T,<j))  =  Xo[n]xl[n  + 

n=l 

1  ^ 

~  ]V  ^  («0®N  +  aix[n  - 

n=l 

{Pox*[n  +  r]  +  /3ix*[n  +  r  - 
=  — ^  x[n]x*  [n  + 

Tl=l 

+  x[n]x*[n  +  T  - 

n=l 

n=\ 

a  3  ^ 

+  Y,  ~  +  T  — 

n=l 

71=1 

+  Y  x[n]x*[n  +  r  - 

n=l 

+  Y  a:[n]a:*[n  +  r  +  ti]e-*(^+“i)" 

n=l— ti 

a;[n]a:*[n  +  T]e-*^ 

n=l—ti 

=  ax{T,4>)  +  bx{T -ti,(j)-wi) 

+  cx{t  +  tu<l)  +  wi)  (35) 
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where 


b  = 
c  = 

This  shows  what  was  hinted  to  in  section  2.2:  that  the  result  of  the  matched  filter  is  the 
sum  of  a  number  of  ambiguity  functions,  scaled  and  shifted  according  to  the  location  and 
velocity  of  the  targets.  However  in  this  case,  a  single  target  has  caused  three  copies  of  the 
ambiguity  function  in  the  result. 

The  complex  phases  on  the  coefficients  a,  b  and  c  are  due  to  the  difierent  time  origins 
of  the  two  signals  which  make  up  the  received  signals.  Components  of  the  signals  are 
delayed  in  time  an  amount  ti,  yet  the  ambiguity  function  for  each  is  calculated  fi:om  the 
same  time  origin.  Like  the  Fourier  transform,  a  time  delay  in  a  signal  causes  a  phase  shift 
in  the  ambiguity  function  of  that  signal. 

The  result  of  squaring  the  matched  filter  output  is: 

+  \4'^\xiT  +  ti,<f>  +  wi)\‘^ 

+  2^{ab*xir,  4>)x*{'r  +  *1,  ^  +  wi)} 

+  2^{ac*x{T,  <I>)x*{t  -tu(f>-wi)} 

+  2?ft{bc*x{'r  +  ti,^  +  wi)x*{T-ti,^-wi)}  (36) 

and 

E{|M(r,<^)|2}  =  E{|ap}E{|x(r,9i)|2} 

+  E{\b\^}E{\x{T-h,cf>-w,)f} 

+  E{|cl2}E{|x(r  +  <i,0  +  «;i)|2} 

+  29?{E{a6*}E{x(r,  +  h,  cj>  +  u;i)}} 

+  25ft{E{ac*}E{x(r,  fl!.)x*(r  -hA- 
+  2SR{E{6c*} 

EWr  +  ti,  ^  +  w{)x*{r  -  ti,  ^  -  tui)}}  (37) 

It  is  assumed  that  the  phase  of  and  /S*,  is  uniformly  distributed  over  the  interval 
[0,27r],  which  implies  that  Efa*}  =  E{j5fc}  =  0  for  A;  =  0  or  1.  Also,  as  the  location  of 
each  of  the  targets  is  independent,  both  aj  and  j3j  are  independent  of  both  and  Pk  for 
3  ^  k. 

These  assumptions  imply  that: 

E{a6*}  =  + 

=  E{/?i}E{|ao|"/3o*}e-'’"^" 

+  EK}E{ai|A|'}e-**^^ 

=  0  (38) 

E{ac*}  =  E{aoail/5o|^e-*(^+“'i)‘i 
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=  E{ao}E{aI|/3op}e-*('^+®i)*i 
=  0 

E{6c*}  =  E{aoat/5o/3ie**i'^-*"’i’-*“i*i} 

=  0 

and  thus: 

E{|M(r,<^)|2}  =  |a|2E{|x(r,<^)|2} 

+  |c|^E{|x(t  +  ti,  ^  +  tni)P} 

Incorporating  information  about  E{|x(r,  ^)p},  it  can  be  seen  that: 


(39) 


(40) 


(41) 


^{\M{T,cf>)\^} 


’  ([Q|2  +  my+|cpj^4 
([6[2  +  Jal!.+|y+|c|^)^4 

^lap+iy+icp^^ 


for  (t,0)  =  (0,0) 
for  (r,^)  =  (ti,wi) 
for  (r,^)  =  (-ii,  -«>i) 
otherwise 


(42) 


This  function  details  the  expected  output  of  the  matched  filter  and  square-law  device 
when  a  single  target  exists  in  the  region  of  interest.  The  function  consists  of  three  spikes, 
at  (0,0),  (ti,u;i)  and  (— ti,— loi),  surrounded  by  a  constant  level  pedestal. 

The  spikes  at  (0,0)  and  (— ti,— Wi)  are  cross-terms  due  to  energy  from  Xo[n]  IpaVing 
into  Xt\p]  and  visa  versa,  and  unless  suppressed,  are  guaranteed  to  produce  false  detections 
when  |A/(t,  ^)P  is  compared  to  a  detection  threshold.  These  cross-terms  are  suppressed 
by  adjusting  the  coefficients  and  /3fc,  whose  magnitudes  are  influenced  heavily  by  the 
beam  pattern  of  the  two  receivers. 

Ideally,  the  magnitudes  of  ai  and  /3o  are  minimised,  so  that  the  only  term  of  significance 
is  h.  Then: 


E{|M(r,<^)|2}  =  |6|"E{|x(r-ii,(^-u;i)|2} 

I  otherwise 

_  f  l«oPl/3iP(l  +  ;^)(r^  (^,0)  =  (ti,ioi) 

\  |aop|/3ip;^cr|  otherwise 

Note  that  this  expression  is  the  same  as  that  for  E{|x('r,  ^)P},  only  shifted  over  (ti,  wi) 
and  scaled  by  |Q!opl/?ip.  Thus  if  the  output  of  the  matched  filter  is  normalised  so  that 
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the  maximum  value  is  tinity,  the  result  is  expected  to  be  a  spike  of  unit  height  at  (ti,  toi) 
surrounded  by  a  constant  level  pedestal  of  value: 


for  large  N  and  (r,<^)  7^  (ti,iwi). 


\ao\W{^  +  i/N)4  ..  1 

lao|21^i|2(l/iV)c74  ~  N 


(44) 


If  the  single  receiver  radar  configuration  of  section  3.1.2  is  used,  the  two  input  signals 
to  the  matched  filter  are  the  same  and  ak  =  ^k-  In  this  case,  it  would  be  impossible  to 
minimise  the  cross-terms  and  the  ability  to  detect  the  target  would  depend  on  the  ability 
to  determine  which  of  the  three  peaks  of  the  matched  filter  corresponded  to  it. 

If  the  received  signal  was  scaled  such  that  |ai|  =  1,  the  result  of  the  matched  filter 
when  a  single  receiver  is  used  is: 


+  I«oPE{|x(t 

+  l«opE{|x(r-fti,<6-hioi)|^}  (45) 


It  is  desired  that  the  peak  at  (ti,  wi)  is  maximised  compared  to  the  peak  at  the  origin, 
as  there  is  concern  that  the  former  peak  may  be  masked  by  the  pedestal  of  the  ambiguity 
function  centred  at  the  origin.  That  is,  the  value  of  |q!o|  is  sought  which  maximises: 


(l  +  |aon 


(46) 


The  required  value  of  |ao|  is  1.  Thus,  when  a  single  receiver  is  used,  it  is  desired  to  fix  the 
attenuation  of  the  antenna  in  the  direction  of  the  transmitter  so  that  the  received  signal 
scattered  firom  the  target  is  equally  as  powerful  as  the  direct  signal. 

Note  that  substituting  |ao|  =  1  into  the  previous  equation  produces  1/2,  thus  when 
a  single  receiver  is  used,  the  maximum  height  of  the  peak  corresponding  to  the  target  is 
3  dB  less  than  the  peak  at  the  origin. 

These  calculations  define  the  expected  value  of  the  matched  filter  output.  An  approx¬ 
imation  of  the  variance  of  this  function  is: 


VAR{lM(r,<^)P}  =  |a|2VAR{lx(r,0)l2} 

+  |6|2  VAR{|x(r-fi,<^-u;i)|2} 

-F|c|2VAR{|x(r  +  ii,<^  +  u;i)|2}  (47) 


where  VAR{|x(t,  ^)|2}  is  given  in  equation  (32). 

When  the  magnitudes  of  ofi  and  are  minimised  and  b  is  the  only  term  of  significance, 
this  expression  simplifies  to: 


VAR{1M(t,<?!.)|2} 


|o:oPlAP(;^  +  ^  for  {T,<f>)  = 

otherwise 


(48) 


Section  4.6  describes  how  these  results  can  be  used  to  determine  the  necessary  threshold 
that  should  be  applied  to  |M(r,  0)l2  in  order  to  detect  targets. 
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4.4  Result  of  the  matched  filter  (general  case) 

If  there  are  K  targets  in  the  region  of  interest,  equations  (33)  and  (34)  can  be  extended 
to  be: 

K 

^o[n]  =  5^  afca;[n  -  (49) 

fc=0 

K 

Xrin]  =  Y.Pkx[n-tk]e-^‘^>‘'^  (50) 

*=0 

where  and  W}-  indicate  the  amount  of  time  delay  and  frequency  shift  produced  by 
the  location  and  velocity  of  the  A:-th  target.  As  the  signal  arriving  directly  from  the 
transmitter  of  opportunity  usually  undergoes  no  delay  in  time  or  shift  in  Doppler  (or  at 
least  all  other  time  delays  and  Doppler  shifts  are  normalised  to  this  signal),  it  is  expected 
that  tk  =  Wk~0  for  one  component  of  these  two  signals  (say,  the  one  for  fc  =  0). 

As  in  the  previous  section,  the  coeflBcients  ak  and  13^  refer  to  the  attenuation  of  the 
component  of  the  signal  reflected  from  the  fc-th  target. 

Applying  equation  (5)  to  the  scenario  described  above,  the  result  of  the  matched  filter 
is: 

1  ^ 

M{T,(f>)  =  —Y^Xo[n]x*[n  +  T]e^^ 

n=l 

=  4  S  S  H  ajl3kx[n  - 

n=lj=0Jfc=0 
x*[n  +  T- 
K  K  N-tj 

^  y  ^  ^  ^  o^j0kx[Ti]x  \n  +  tj  tfc  "I" 

j—0  k=0n=l—tj 

Q—i(wjn+Wjtj—wi,n—wictj—WjT+(lm+<j>tj) 

^  K  K 

=  aj0ke~^^  (^+™j  -'U)k)tj-wkT) 

^  j=o k=o 
N-tj 

x[n]x*[n  +  (r  +  tj  — 

n=l—tj 
K  K 

=  +  tj-tk,<i>3-Wj-Wk)  (51) 

j=0k=0 

where  jj^k  = 

The  result  of  squaring  the  matched  filter  output  is: 

=  42E  E  E  E7fex,ife.7fc3,/b4 

A:i=0  A:2=0  k3=0  fc4=0 
X(x  +  tki  -tk2,(l>  +  Wkj  -  Wk^) 

X* {t  +  43  -tki,(f>  +  Wks  -  Wk^ )  (52) 
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and 


K  K  K  K 

E{|M{r,^)P}  =  E  E  EE{7«.,/i,7U> 

A:i=0  fc2=0  fc3=0  fc4=0 

E  {x(t  +  4i  -tkiA  + 

X*i'r  +  tks  -tk4,<t>  +  Wks  -  1«A:4)} 


(53) 


where 


E{'yfci,fc2Tfc3,fc4}  ^i-^kiPk^^ks^ki} 

g-t((</>+iu*!i  -Wk2  )tki  )tk3  -(®*2 

Again  it  is  assumed  that  E{q!*}  =  E{/3fc}  =  0  and  that  both  aj  and  j3j  are  independent 
of  both  ak  and  Pk  for  j  ^  k. 

These  assumptions  imply  that: 


E{TA:i  ,fc2  Ifka  ,k4  }  ^ki—ka  ^k2—k4  E{Q:fci  Pk2  ^*3  ^*4  } 

Q-Ki^+Wki  -Wk2  )tki  -(.<I>+Wk2 -Wk^  )tk3  ) 

where  Sj  is  1  for  j  =  0  and  0  otherwise. 

Equation  (53)  can  thus  be  written  as: 

K  K 
j=0fc=0 

E{|x(t  +  tj  -tk,(l>  +  Wj  -  u;fc)p} 

=  ^EEEdoiflAp} 

j=0  k=0 

E{|x('r  +  tj  -tk,<f>  +  Wj  -  Wk)\^} 

This  equation  shows  that  the  square  of  the  matched  filter  is  the  sum  of  up  to  {K  + 1)^ 
copies  of  the  ambiguity  function,  which  appear  as  (if  +  1)^  peaks  across  the  (r,  4>)  plane. 
Of  these,  only  K  correspond  to  actual  targets  and  the  remainder  are  the  undesired  cross¬ 
terms. 

As  for  the  single  target  case,  these  cross-terms  are  reduced  by  adjusting  the  magnitudes 
of  ak  and  Pk-  Note  that  effective  reduction  in  these  cross-terms  is  possible  by  manipulating 
only  one  of  these  sets  of  coefficients.  For  example,  if  the  receiver  collecting  the  direct  signal 
substantially  rejects  signals  from  any  other  direction,  |Q!fc|  will  be  approximately  zero  for 
all  k  but  one.  If  it  is  |ao|  that  is  non-zero,  and  the  direct  signal  is  indeed  direct  (that  is, 
to  =  =  0),  then  equation  (53)  becomes: 

E{|M{t,,^)P}  =  Ey;E{|aoP|ftP} 

/b=0 

^{\x{t  -tk,<p-Wk)\'^}  (56) 


(54) 


(55) 
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which  corresponds  to  a  peak  at  the  origin  and  K  other  peaks,  each  located  at  a  time  delay, 
tk,  and  Doppler  shift,  Wk,  corresponding  to  the  fe-th  target. 

In  the  configuration  using  a  single  antenna,  =  Pk  and  Xo\p\  =  a;r[Ti].  In  this  case 
equation  (53)  becomes: 

K  K 

E{|M(r,^)P}  =  j4EEB{l“iPKI"} 

j=0k=0 

H\x{r  +  tj-tk,4>  +  Wj-Wk)\^}  (57) 

Thus  for  the  single  receiver  configuration,  it  is  much  more  difficult  to  reduce  the  cross¬ 
terms  in  the  matched  filter  output. 


4.5  Resolution  issues  re-visited 

Section  2.3  provided  a  general  discussion  on  how  the  shape  of  the  radar  waveform,  and 
more  specifically  the  shape  of  its  ambiguity  function,  can  affect  the  ability  of  the  radar 
to  resolve  target  location  and  velocity.  This  section  examines  the  results  of  the  previous 
three  sections  to  draw  more  specific  conclusions  on  the  performance  of  the  proposed  radar 
system  and  waveforms  of  opportunity. 

It  is  evident  that  the  output  of  the  matched  filter  is  the  combination  of  a  number  of 
ambiguity  functions,  each  shifted  to  the  time  delay  and  Doppler  shift  corresponding  to 
a  target  (or  false  target  in  the  case  of  cross-terms),  and  scaled  by  some  amount.  Each 
ambiguity  function  consists  of  a  peak  at  the  origin  smrounded  by  a  constant  level  pedestal. 
This  result  suggests  two  possible  problems: 

•  That  the  peaks  due  to  cross-terms  produce  false  target  detections,  and 

•  That  a  copy  of  the  ambiguity  function  due  to  one  target  is  scaled  such  that  the 
constant  level  pedestal  away  from  its  origin  is  larger  than  the  peak  of  the  ambiguity 
frmction  due  to  another  target,  rendering  the  second  tajget  undetectable. 

As  mentioned  already,  the  solution  to  the  first  problem  is  to  adjust  the  antenna  gain 
pattern.  Note  that  there  need  not  be  perfect  attenuation  of  signals  arriving  from  directions 
other  that  the  direction  of  interest  for  the  cross-terms  to  be  avoided.  For  example,  if  there 
exists  two  targets  and  the  powers  in  the  received  signals  are  normalised  so  that  the  direct 
wave  has  a  power  of  unity,  the  scattered  signals  will  typically  contain  power  much  less 
than  imity,  that  is,  |ao|  in  equation  (49)  is  1  and  |/3i|  and  \^\  in  equation  (50)  are  less 
than  1.  The  peaks  in  the  matched  filter  output  corresponding  to  the  two  targets  will  have 
height  |/3ip  and  |/%p  while  the  cross-terms  will  have  height  |/0ip  x  1/^2 1^,  which  is  far  less. 

Although  cross-terms  produced  by  the  interaction  of  retmrns  from  multiple  targets  are 
often  negligible,  the  interaction  between  the  amount  of  the  direct  signal  in  both  Xr[n] 
and  rcojn]  usually  produces  a  peak  in  the  origin  of  the  matched  filter  output  of  significant 
magnitude.  This  cross-term  is  usually  the  cause  of  the  second  problem  of  the  two  listed 
above. 
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The  power  of  the  direct  signal  is: 


p 

°  (4<r)V§ 

where  Pt  is  the  power  used  by  the  transmitter,  A  is  the  wavelength  of  the  carrier  frequency 
and  To  is  the  distance  from  the  transmitter  to  the  receiver.  The  power  of  a  signal  scattered 
from  a  target  and  received  at  the  radar  is; 


Pr  = 


Pt>?o 

(47r)3(rir2)^ 


(59) 


where  r\  is  the  distance  from  the  transmitter  to  the  target  and  r2  is  the  distance  from  the 
target  to  the  receiver. 

Typically,  P^  is  significantly  less  than  Pg  so  even  with  substantial  attenuation  of  the 
direct  signal,  the  coefficient  |^|  in  equation  (34)  is  significant  and  the  matched  filter 
output  contains  the  peak  at  the  origin  that  was  mentioned  earlier. 

The  ratio  of  the  power  of  the  two  received  signals  is: 

Po  47r(rir2)2 

This  value  is  the  difference  in  height  of  the  two  peaks  in  the  matched  filter  output  that 
occur  if  the  antenna  at  the  receiver  was  onmi-direction.  If  this  value  is  more  negative  than 
the  level  of  the  pedestal  of  the  ambiguity  function  centred  at  the  origin,  the  target  is  lost. 
Attenuation  towards  the  transmitter  reduces  Pg  and  increases  the  height  of  the  peak  due 
to  the  target  compared  to  the  peak  at  the  origin.  The  amount  of  attenuation  required 
is  the  amoimt  which  brings  the  height  of  the  peak  due  to  the  target  above  the  pedestal 
around  the  peak  at  the  origin. 

The  constant  level  pedestal  surrounding  the  various  peaks  in  the  matched  filter  output 
is  often  called  the  clutter. 


4.6  Target  detection 

As  mentioned  in  section  2.1,  optimally  detecting  targets  involves  comparing  the  output 
of  the  matched  filter,  |M(t,  ^)P,  against  a  single  threshold  and  determining  at  which  values 
of  T  and  (j)  the  threshold  is  exceeded.  An  obvious  question  is  what  to  set  the  threshold  to. 

The  previous  sections  of  this  chapter  showed  that  the  output  of  the  matched  filter 
and  square-law  device  is  a  number  of  peaJcs  corresponding  to  the  returns  from  targets 
surrounded  by  a  constant  level  pedestal.  Thus  the  detection  threshold  is  required  to 
separate  the  peaks  from  the  pedestal. 

Detection  theory  provides  a  number  of  methods  for  which  to  optimally  determine 
this  threshold,  however  the  choices  are  drastically  reduced  for  the  radar  problem  as  the 
likelihood  of  a  target  being  present  at  a  given  time  delay  or  Doppler  shift  cannot  be  known. 
One  method  which  is  often  used  is  the  Neyman-Pearson  detector  [27]. 
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The  Neyman-Pearson  detector  involves  calculating  the  detection  threshold  so  that  the 
probability  of  falsely  determining  the  presence  of  a  target  is  fixed.  For  the  receiver  output 
described  above  this  involves  calculating  the  threshold  so  that  a  fixed  number  of  (r,  <f)) 
pairs  which  occur  in  the  pedestal  produce  values  of  |M(r,  which  exceed  the  threshold. 

This  calculation  requires  the  probabilistic  distribution  of  the  values  of  the  output  of 
the  matched  filter  receiver  which  lie  in  the  constant  level  pedestal.  This  could  be  estimated 
by  sampling  the  matched  filter  output  at  values  of  r  and  <f)  which  are  known  not  to  contain 
any  peaks  due  to  targets. 

Alternatively,  the  distribution  could  be  assumed  to  be  a  particular  type,  such  as  Gaus¬ 
sian  or  chi-squared,  which  is  totally  described  by  the  two  moments  (mean  and  variance) 
calculated  in  the  previous  sections.  For  example,  if  the  distribution  is  assumed  Gaussian 
and  the  probability  of  false  alarm  is  set  to  1  in  100,  the  detection  threshold  would  be  3 
standard  deviations  above  the  mean.  If  there  existed  a  single  target  in  the  region  of  inter¬ 
est  and  the  antennas  were  such  that  the  amount  of  direct  signal  contributing  to  Xr[n]  and 
the  amount  of  scattered  signal  contributing  to  Xo\t]  was  negligible,  the  mean  of  the  output 
of  the  matched  filter  receiver  is  given  in  equation  (43)  and  its  variance  in  equation  (48). 
If  the  output  is  normalised  so  that  the  peak  at  [ti ,  toi )  is  unity,  the  necessary  detection 
threshold  would  be: 


Often  this  detection  step  is  not  the  final  processing  performed  on  the  received  signals. 
Tracking  and  the  algorithms  mentioned  in  section  3.2  which  consider  how  target  detections 
vary  with  time  are  also  used  to  minimise  false  detections.  For  this  reason,  the  probability 
of  false  detection  used  when  calculating  the  above  threshold  is  often  fixed  to  be  slightly 
more  than  what  is  ultimately  desired. 
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Chapter  5 


Simulations  and  Examples 


The  previous  chapters  provide  insight  as  to  how  it  may  be  possible  to  use  transmissions 
such  as  radio  and  television  broadcasts  as  the  waveform  for  passive  radar.  This  chapter 
provides  a  number  of  simulations  of  this  concept  and  discussions  on  how  the  ideas  presented 
in  this  report  could  be  used. 

These  simulations  provide  some  confirmation  of  the  results  attained  in  chapter  4  and 
an  initial  picture  of  the  performance  of  the  suggested  passive  radar. 


5.1  FM  radio 

One  possibility  for  the  proposed  passive  radar  is  to  inconspicuously  monitor  ships 
and  aircraft  arriving  at  and  departing  from  a  particular  city  or  location  of  interest.  The 
example  city  chosen  for  the  first  simulation  was  Darwin,  the  largest  port  along  Australia’s 
northern  coast  and  the  transmitter  of  opportunity  was  a  local  PM  radio  station  located 
on  the  city’s  edge. 

The  receiver  was  located  62  km  away.  It  is  worth  noting  again  that  the  radar  emits  no 
signal  energy  and  could  therefore  be  hidden. 

This  simulation  examines  the  ability  to  detect  a  single  ship  which  is  approaching  Dar¬ 
win.  An  illustration  of  the  scenario  is  shown  in  figme  5.1  and  the  parameters  of  simulation 
are  listed  in  the  following  table: 


Simulated  scenario: 

Distance  from  Transmitter  to  Receiver  (ro) 

62.0  km 

Distance  from  Transmitter  to  Target  (ri) 

21.4  km 

Distance  from  Target  to  Receiver  (r2) 

73.5  km 

Target  velocity  (v) 

48  km/hr 

Target  RCS  (a) 

250  m2 

Transmitter  carrier  frequency  (/c) 

103.3  MHz 

Receiver  characteristics: 

Attenuation  directly  towards  the  transmitter  {^q) 

-50  dB 

Radar  integration  time  (T) 

2.8  s 

Receiver  bandwidth  (B) 

4  kHz 
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It  is  assumed  that  a  linear  array  of  antennas  was  being  used  to  collect  the  incoming 
waveforms.  In  theory  such  arrays  are  capable  of  totally  rejecting  signals  arriving  from  a 
particular  direction  however  in  practice  the  depth  of  any  null  is  governed  by  how  well  the 
array  is  calibrated.  As  emphasised  already  in  earlier  chapters,  the  ability  to  attenuate  the 
direct  signal  is  very  important  to  the  performance  of  the  radar. 

As  with  all  the  simulations  in  this  chapter,  the  received  signals  were  considered  digital, 
having  been  sampled  upon  being  received.  Thus,  the  result  of  the  matched  filter  was 
evaluated  at  intervals  of  1/5  along  the  t  axis  and  1/T  along  the  (f>  axis  (this  process  is 
explained  further  in  section  2.5).  This  implies  the  length  of  a  resulting  range-Doppler  cell 
was  c/4000  =  75  km  along  the  axis  corresponding  to  total  range  and  cf{Tfc)  =  3.73  km/hr 
along  the  axis  corresponding  to  velocity.  Thus  the  parameters  of  the  radar  defined  here 
have  the  potential  to  provide  satisfactory  resolution  in  target  velocity  but  not  in  range  for 
the  given  scenario.  This  is  illustrated  in  figure  5.1  which  shows  ellipses  marking  the  edge 
of  the  first  and  second  range  cells. 

The  ambiguity  function  of  the  waveform  of  opportunity  which  was  used  in  this  simu¬ 
lation  is  shown  in  figure  5.2.  This  waveform  was  an  actual  recording  from  an  FM  radio 
transmission. 

The  radar  waveform  was  delayed  in  time  and  shifted  in  Doppler  the  appropriate  amount 
to  simulate  the  location  and  velocity  of  the  target.  This  signal  and  the  original  waveform 
were  both  scaled  to  simulate  the  path  loss  due  to  distance  travelled  and  combined  to  form 
the  two  received  signals  Xo[n\  and  Xr[n\.  These  signals  became  the  input  to  the  matched 
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Figure  5.1:  First  proposed  simulation.  The  dotted  ellipses  mark  the  extent  of  the  first  and 
second  range  cells. 
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filter,  as  dictated  in  equation  (5).  The  output  of  this  filter  and  the  result  of  this  first 
simulation  is  shown  in  figme  5.3. 

In  this  figure  the  peak  due  to  the  target  can  be  just  seen  above  the  clutter  floor  caused 
by  the  direct  signal  and  the  copy  of  the  ambiguity  function  located  at  the  origin.  The  level 
of  this  clutter  floor  is  expected  to  be  1/(T5)  =  8.93  x  10“®  =  —40.5  dB  and  the  ratio  of 
power  in  the  scattered  signal  to  that  in  the  direct  signal  is: 

Pr  1  250  X  620002 

Po  ~  47r  (21400  X  73500)2 
=  3.09x10"® 

=  -75.1  dB 

Prom  this  it  can  be  seen  that  without  the  attenuation  of  the  array  in  the  direction  of 
the  transmitter,  the  power  in  the  scattered  signal  would  be  dominated  by  the  clutter 
floor.  With  the  attenuation  in  place  however,  the  difference  between  the  powers  becomes 
—75.1  +  50  =  —25.1  dB,  pushing  the  height  of  the  peak  due  to  the  target  above  the  clutter 
floor.  These  results  can  be  seen  in  figure  5.3. 

Using  equation  (48)  and  the  fact  that  the  matched  filter  output  has  been  normalised, 
the  variance  of  the  clutter  floor  is: 

_ ? _ I _ ^ 2  28  X  in~® 

(2.8  X  4000)2  ^  (2.8  x  4000)^ 


Figure  5.2:  Ambiguity  function  of  the  waveform  of  opportunity  for  the  first  simulation. 
The  waveform  had  length  T  =  2.8  s  and  bandwidth  B  =  A  kHz. 
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Figure  5.3:  Result  of  the  matched  filter  for  the  first  simulation. 
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Figure  5-4:  Result  of  the  matched  filter  after  thresholding  at  3  standard  deviations  above 
the  expected  clutter  floor  (first  simulation). 


Thus,  three  standard  deviations  above  the  clutter  floor  comes  to  -32.7  dB.  Thresholding 
the  output  of  the  matched  filter  at  this  value  is  illustrated  in  figure  5.4. 

This  figure  clearly  shows  both  the  peak  due  to  the  direct  signal  and  the  peak  due  to  the 
target.  Also  apparent  in  this  figure  are  six  range-Doppler  cells  where  the  matched  filter 
output  exceeded  the  threshold  at  values  of  r  and  (p  for  which  there  was  no  target.  The 
figure  contains  640  range-Doppler  cells  corresponding  to  the  clutter  floor,  so  the  probability 
of  false  alarm  appears  to  be  6/640  =  0.009,  close  to  the  expected  1  percent. 

Prom  these  results  it  seems  the  suggested  configuration  was  able  to  successfully  de¬ 
tect  the  target  and  calculate  its  total  range  and  rate  of  change  of  total  range^  although 
range  resolution  was  particularly  poor.  In  an  effort  to  improve  upon  range  resolution,  the 
previous  simulation  was  executed  a  second  time  using  a  signal  of  bandwidth  B  =  48  kHz. 

In  this  second  simulation,  all  parameters  remained  the  same  except  the  bandwidth  of 
the  receivers.  The  length  of  a  range-Doppler  cell  along  the  range  axis  for  this  simulation 
IS  c/B  =  6.25  km,  thus  it  is  expected  to  be  able  to  resolve  the  location  of  targets  on  a 
much  finer  grid.  This  is  illustrated  in  figure  5.5. 

The  ambiguity  function  of  the  waveform  used  in  this  second  simulation  is  shown  in 
figure  5.6.  It  should  be  noted  that  this  signal  was  simulated  by  frequency  modulating  a 
sample  of  music  as  a  receiver  with  a  bandwidth  of  up  to  48  kHz  was  unavailable  during  the 
course  of  preparing  this  report.  Nevertheless,  the  ambiguity  function  appears  typical  of 

‘‘Conversion  from  total  range  and  rate  of  change  of  total  range  to  target  location  and  velocity  is  provided 
in  section  3.2. 
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Degrees  east 


Figure  5.5:  Second  proposed  simulation.  The  dotted  ellipses  mark  the  extent  of  the  ranges 
cells. 


that  generated  from  a  transmission  of  opportunity  and  contains  the  expected  constant  level 
pedestal  of  ll{TB)  =  -51.3  dB  {more  than  10  dB  lower  than  the  pedestal  of  figure  5.1). 

The  result  of  the  matched  filter  can  be  seen  in  figure  5.7.  As  the  location  of  the 
target  and  amount  of  attenuation  of  the  receiver  towards  the  transmitter  has  not  changed, 
the  difference  of  heights  of  the  two  peaks  in  \M{t,4>)\^  is  still  -25.1  dB.  The  increase  in 
bandwidth  however,  has  caused  the  clutter  fioor  around  the  peaks  to  drop  10.8  dB.  This 
implies  that  the  attenuation  of  the  receivers  in  the  direction  of  the  direct  signal  could  have 
been  10.8  dB  worse  and  the  target  would  still  have  been  as  visible  as  it  was  in  the  first 
simulation. 

This  figure  also  shows  that  the  location  of  the  target  causes  a  peak  in  the  fifth  range 
cell,  as  opposed  to  the  first,  which  was  the  case  in  the  first  simulation.  Thus  it  is  possible 
to  resolve  to  better  accTiracy  the  location  of  the  target. 

For  this  second  simulation,  three  standard  deviations  above  the  clutter  floor  corre¬ 
sponds  to  —44.1  dB.  The  result  of  the  matched  filter  thresholded  at  this  value  is  shown 
in  figure  5.8.  As  can  be  seen  from  this  figure,  the  peak  corresponding  to  the  target  sits 
strong  above  the  threshold,  however  there  appears  to  be  a  much  greater  amount  of  the 
|M(r,  surface  which  also  sits  above  it.  Much  of  these  detections  are  due  to  the  wide 
peak  at  the  origin.  The  transmissions  of  opportunity  are  not  exactly  white  noise  and  the 
peak  at  the  origin  of  their  ambiguity  function  has  finite  width. 

One  final  note  on  the  benefits  of  increasing  receiver  bandwidth  is  shown  in  figure  5.9. 
Part  (a)  of  this  figure  shows,  for  each  location  on  the  map,  the  amount  of  energy  arriving 
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Figure  5.6:  Ambiguity  function  of  the  waveform  of  opportunity  for  the  second  simulation. 
The  waveform  had  length  T  =  2.8  s  and  bandwidth  J5  =  48  kHz. 

at  the  receiver  due  to  the  scattering  off  a  target  at  that  location.  Part  (b)  shows  the 
same  information  as  part  (a),  but  only  for  those  locations  which  produce  signals  of  energy 
—44.1  —  50  =  —94.1  dB  and  —32.7  —  50  =  —82.7  dB.  That  is,  the  two  ellipses  of  part  (b) 
mark  the  edge  of  detection  for  the  receivers  of  the  two  simulations  (the  inside  ellipse 
corresponds  to  the  detection  threshold  used  in  the  first  simulation).  Any  target  located 
beyond  an  ellipse  produces  a  scattered  signal  whose  energy  would  not  be  enough  for  the 
peak  in  the  matched  filter  output  to  appear  above  the  given  detection  threshold. 

As  can  be  seen  from  this  figure,  the  receiver  using  the  higher  bandwidth  can  detect 
targets  at  a  much  further  range. 


5.2  Television  broadcasts 

The  previous  section  showed  that  by  increasing  the  bandwidth  of  the  receiver  to  48  kHz, 
the  radar  was  able  to  resolve  the  location  of  the  target  to  within  approximately  6  km.  In 
order  to  resolve  targets  to  a  finer  resolution,  the  required  bandwidth  is  greater.  Of  the 
three  t3rpes  of  transmissions  of  opportunity  investigated  in  this  report,  the  only  ones  which 
contain  energy  distributed  across  such  large  bandwidths  are  television  broadcasts. 

As  mentioned  in  section  3.3,  television  broadcasts  support  a  bandwidth  of  6  MHz.  If 
the  bandwidth  of  the  receiver  was  half  of  this,  the  location  of  targets  could  be  determined 
to  within  c/(3  x  10®)  =  100  m.  Television  broadcasts  however,  contain  a  dominant  pulse 
occurring  every  64  fis  which  produce  ridges  64  )us  apart  in  the  ambiguity  function  and 
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Figure  5. 1:  Result  of  the  matched  filter  for  the  second  simulation. 
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Figure  5.8:  Result  of  the  matched  filter  after  thresholding  at  3  standard  deviations  above 
the  expected  clutter  floor  (second  simulation). 


multiple  peaks  in  the  matched  filter  output  when  a  single  target  is  present  in  the  region 
of  interest. 

A  time  delay  of  64  //s  corresponds  to  a  total  range  of  c  x  64  x  10“®  =  19.2  km.  Thus  a 
taxget  whose  total  range  was  19.2  +  d  km  could  be  mistaken  for  a  target  whose  total  range 
was  only  d  km.  One  solution  to  this  ambiguity  is  to  scan  only  for  taxgets  of  total  range 
less  than  19.2  +  ro  km  and  to  select  a  time-bandwidth  product  such  that  the  returns  from 
targets  beyond  this  range  fall  below  the  detection  threshold. 

The  scenario  chosen  for  the  simulation  utilising  television  broadcasts  consists  of  the 
receiver  located  r2  =  2  km  from  an  airport  of  interest  and  tq  =  50  km  from  a  television 
broadcast  tower.  The  relative  locations  of  the  receiver,  transmitter  and  airport  is  shown 
in  figure  5.10.  Also  shown  on  this  figure  is,  for  each  (a;,  y)  location,  the  amount  of  energy 
arriving  at  the  receiver  due  to  the  scattering  off  a  target  at  that  location,  and  an  ellipse 
corresponding  to  a  total  bistatic  range  of  50  -1- 19.2  =  69.2  km.  The  RCS  of  the  target  was 
assumed  to  be  150  m^. 

In  order  to  minimise  the  detection  of  targets  beyond  the  ambiguous  range,  the  receiver 
need  only  be  accurate  enough  to  detect  scattered  signals  with  power  greater  than  approx¬ 
imately  -70  dB  lower  than  that  from  the  direct  signal.  This  allows  the  receiver  to  scan  for 
targets  to,  at  best,  approximately  10  km  away  from  the  airport. 

If  the  receiver  antenna  provided  20  dB  of  attenuation  towards  the  transmitter,  the 
required  detection  threshold  would  be  about  -50  dB.  Given  that  the  receiver  uses  a  band¬ 
width  of  3  MHz,  the  signal  duration  required  to  move  the  clutter  fioor  below  this  threshold 
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Degrees  east 


Figure  5.9:  (a)  Power  received  from  the  scattering  from  a  target  given  its  location  and  (b) 
threshold  for  three  standard  deviations  above  the  noise  floor  for  the  first  two  simulations, 
given  50  dB  attenuation  towards  the  transmitter. 
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km  east 


Figure  5.10:  Power  received  from  the  scattering  from  a  target  given  its  location  for  the 
third  simulation.  The  dashed  ellipse  marks  a  total  bistatic  range  of  69.2  km. 

is  approximately  T  =  x  3  x  10®)  =  0.03  s. 

Given  this  integration  time,  the  velocity  resolution  (length  of  a  cell  in  the  matched 
filter  output  along  the  Doppler  axis)  is  45  m/s=162  km/hr. 

A  more  desirable  solution  to  the  problem  of  ambiguous  range  is  to  analyse  the  location 
and  velocity  of  targets  over  time,  as  mentioned  in  section  3.2.  Targets  should  vary  their 
range  over  successive  integration  intervals  as  dictated  by  their  velocity,  and  as  integration 
time  is  small,  the  most  likely  of  the  possible  ranges  of  a  particular  target  could  be  estimated 
in  a  comparatively  short  period. 


5.3  Measuring  the  distance  to  the  moon 

The  simulations  and  examples  provided  thus  far  have  shown  how  the  proposed  passive 
radar  can  be  used  to  locate  aircraft  and  ships  moving  relatively  close  to  the  receiver.  The 
same  techniques  however,  could  be  used  to  undertake  a  more  novel  project:  estimating 
the  distance  to  the  moon.  For  this  task,  the  transmitter  of  opportunity  would  be  the  sun. 

To  measure  the  distance  to  the  moon,  both  the  moon  and  the  sun  are  required  to  be 
in  the  visible  sky,  so  that  both  the  direct  signal  and  signal  scattered  fi:om  the  moon  are 
available  to  the  receiver.  An  example  of  this  is  during  the  early  evening  when  the  moon 
is  near  full.  The  basic  geometry  of  this  experiment  is  illustrated  in  figure  5.11. 

Most  of  the  energy  transmitted  by  the  sun  occurs  at  the  firequencies  corresponding  to 
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Figure  5.11:  Proposed  configuration  for  measuring  the  range  and  velocity  of  the  moon. 


visible  light,  although  there  is  still  much  energy  which  occurs  at  lower  frequencies,  such  as 
those  at  the  high  end  of  the  radio  spectrum.  The  U.S.  National  Telecommunications  and 
Information  Administration  (NTIA)  allocates  frequencies  between  250  and  252  GHz  for 
passively  monitoring  radiation  from  space,  so  a  carrier  frequency  within  this  band  seems  a 
reasonable  choice.  It  is  expected  that  the  transmissions  from  the  sun  at  these  frequencies 
are  very  similar  to  white  noise  over  realistic  receiver  bandwidths. 

It  is  assumed  that  the  distance  from  the  earth  to  the  sun  is  far  greater  than  that  from 
the  earth  to  the  moon,  so  n  «  tq  =  149.6  x  10®  m.  This  implies  that  the  distance  to 
the  moon,  r2,  can  be  simply  calculated  using  ctd  (where  ta  is  the  time  delay  of  the  peak 
in  the  matched  filter  corresponding  to  the  return  from  the  moon),  and  that  an  accurate 
measure  of  the  angle  between  the  direction  of  arrival  of  the  two  received  signals,  9,  is  not 
important. 

The  distance  to  the  moon  is  384.4  x  10®  m,  and  a  very  approximate  cross-section  of 
the  moon  is  7ri?^  =  tt  x  (3.476  x  10®)®  =  38  x  10^®  m®  (where  Rm  is  the  radius  of  the 
moon).  Thus  the  ratio  of  power  in  the  scattered  signal  to  that  in  the  direct  signal  is: 

Pr  _  1  38  X  10^®  X  (149.6  X  10®)® 

Po  ~  47r  (149.6  X  10®  x  384.4  x  10®)® 

=  20.4  X  10“® 

=  -47  dB 

It  is  assumed  that  the  moon  is  the  only  dominant  reflector  in  the  sky,  so  that  if  the 
single  receiver  configuration  of  section  3.1.2  was  used,  the  only  cross-terms  would  occur  at 
(0, 0)  and  {—t^,  —Wm),  where  tm  and  Wm  are  the  time  time  delay  and  Doppler  shift  caused 
by  the  moon.  Section  4.3  explains  that  for  this  configuration,  the  peak  corresponding  to 
the  target  is  maximised  when  the  attenuation  towards  the  transmitter  is  such  that  the 
strength  of  the  signal  scattered  from  the  target  is  equal  to  that  of  the  direct  signal. 

Thus,  it  is  desired  that  the  receiver  attenuate  signals  arriving  directly  from  the  sun  by 
—47  dB.  If  this  occurs,  the  peaJs  in  the  output  of  the  matched  filter  at  {tm,Wm)  is  3  dB 
lower  than  the  level  of  the  peak  at  the  origin  and  thus  sits  above  the  clutter  floor  for  any 
reasonable  time-bandwidth  product. 

As  the  amount  of  attenuation  towards  the  sun  varies  from  —47  dB,  the  level  of  the  peak 
in  the  matched  filter  corresponding  to  the  moon  decreases,  as  dictated  by  equation  (46). 
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For  example,  if  the  antenna  was  able  to  attenuate  the  direct  signal  by  only  30  dB,  the 
height  of  the  peak  corresponding  to  the  moon  would  be: 


|10(47-3O)/lO|2 
1  +  |10(47-30)/10|4 


=  3.98  X  10-^  =  -34  dB 


The  time-bandwidth  product  of  the  radar  must  be  enough  to  produce  a  clutter  floor  whose 
level  is  lower  than  this.  If  a  receiver  with  B  =  500  kHz  bandwidth  was  used  and  a  signal  of 
length  T  =  4  s  was  recorded,  the  clutter  floor  due  to  the  central  peak  would  be  at  —63  dB 
and  the  threshold  of  three  standard  deviations  above  this  would  correspond  to  —56  dB. 
Thus  the  peak  revealing  the  time  delay  and  Doppler  shift  produced  by  the  moon  would 
be  easily  detected. 

The  distance  to  the  moon  can  be  measured  using  a  single  omni-directional  antenna. 
In  this  case,  the  peak  at  {tm,Wm)  would  be  at  —97  dB.  Even  a  small  bandwidth  receiver 
can  produce  a  clutter  pedestal  lower  than  this  if  the  integration  time  is  sufficient. 
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Chapter  6 

Conclusions 


The  intention  of  this  report  was  to  investigate  the  suitability  of  a  number  of  broadcast 
band  transmitters  of  opportunity  for  a  passive  radar.  It  is  believed  that  a  significant 
amount  of  work  has  been  completed  in  this  task. 

This  report  has  shown  that  signals  typical  of  AM  radio,  FM  radio  and  television  broad¬ 
casting  are  indeed  suitable  radar  waveforms  provided  that  a  number  of  other  parameters 
of  the  radar  are  adequately  tuned.  These  parameters  are: 

•  The  time-bandwidth  product,  TB  (where  T  is  the  radar  integration  time  and  B  is 
the  receiver  bandwidth),  and 

•  The  gain  pattern  of  the  receiver  antenna 


Of  these,  the  receiver  bandwidth  is  severely  limited  by  the  choice  of  transmitter  of  oppor¬ 
tunity. 

The  most  significant  problem  of  using  the  proposed  transmissions  for  passive  radar  is 
the  possibility  of  the  clutter  floor  due  to  the  signal  received  directly  from  the  transmitter 
masking  the  signal  scattered  fi:om  a  target.  This  is  controlled  by  ensuring  the  time- 
bandwidth  product  is  large  and  by  attenuating  signals  arriving  firom  the  direction  of  the 
transmitter. 

Ideally,  two  receivers  (and  corresponding  directive  antenna)  are  required  to  detect 
targets,  one  collecting  the  direct  signal  and  one  scanning  for  targets  in  a  bearing  of  interest. 
It  has  been  shown  however,  that  it  is  possible  for  the  passive  radar  to  operate  using  a  single 
receiver  (and  antenna),  although  the  performance  is  limited  and  more  restrictions  exist  on 
the  three  parameters  mentioned  earlier. 

Of  the  three  t3q)es  of  transmissions  of  opportunity  examined,  those  firom  FM  radio 
appeared  the  most  suitable  as  radar  waveforms.  These  provided  a  signal  closest  to  white- 
noise  and  up  to  a  bandwidth  of  75  kHz  (a  range  resolution  of  up  to  4  km  cells). 

Although  this  report  provides  a  substantial  investigation  into  the  suitability  of  wave¬ 
forms  of  opportunity,  there  exists  many  avenues  of  further  study  in  the  field.  Further  work 
includes: 
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•  Including  noise  and  signals  scattered  from  the  environment  in  the  model  of  the 
received  signal.  In  this  report,  it  was  assumed  that  these  effects  were  negligible 
compared  to  the  clutter  due  the  direct  transmitted  signal. 

•  An  investigation  into  the  use  of  alternate  detection  strategies.  This  report  has  con¬ 
sidered  only  applying  a  single  threshold  to  the  result  of  a  matched  filter,  however 
more  ad  hoc  algorithms  may  prove  more  successful  given  the  particular  waveforms. 

•  An  investigation  into  the  improvement  gained  by  processing  sequences  of  successive 
detections.  This  relates  to  tracking  and  is  considered  beyond  the  scope  of  this  report. 

•  Analysing  the  transmissions  of  opportunity  other  than  the  three  types  proposed  in 
this  report.  Many  more  surely  exist,  including  mobile  communications  and  satel¬ 
lite  transmissions,  which  may  provide  more  bandwidth  and  better  self-correlation 
characteristics. 

The  most  obvious  continuation  of  the  work  presented  in  this  report,  however,  is  to 
physically  implement  the  proposed  passive  radar. 
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Appendix  A:  Detail  of  Mathematical  Analysis 


A.l  Resolving  range  in  bistatic  radar 

Figure  3.1  illustrates  a  typical  bistatic  radar.  The  direction  of  the  main  lobe  of  the 
reciever  provides  9  and  the  matched  filter  provides  the  time  delay  between  the  direct  and 
the  reflected  signal,  td,  such  that; 

ri  +  r2  -  ro  =  ctd  (Al) 

where  c  is  the  speed  of  light.  This  implies: 

Ti  =  {ctd  +  »'o  -  ’”2)^ 

=  (?td  +  »'o  +  ^2  +  2ctdro  -  2ctdr2  -  2ror2  (A2) 


Prom  the  geometry  of  the  triangle  ro  —  ri  —  r2  [28]: 


=  Tq  +  rl  —  2ror2  cos  9 


Combining  equations  (A2)  and  (A3): 

cHd  +  2ctdro  -  2ctdr2  -  2ror2  =  —2ror2  cos  9 
=>  2{ctd  +  ro  -  ro  cos  9)r2  =  cHl  +  2ctdro 

_  ctdjctd  +  2ro) 

2  {ctd  +  ro(l  —  cos  9)) 


A. 2  Resolving  velocity  in  bistatic  radar 

As  in  the  previous  section,  the  geometry  assumed  in  calculating  the  velocity  of  a  target 
is  shown  in  figure  3.1. 

The  matched  filter  provides  the  amount  of  firequency  shift,  Wd,  between  the  two  received 
signals,  which  relates  to  the  rate  of  change  of  the  total  range  of  the  target,  ri  +  r2,  in  the 
following  manner: 

dt  ^  dt  2nfc 


where  fc  is  the  carrier  frequency  of  the  radar  transmission. 

Differentiating  equation  (A3)  with  respect  to  t  provides: 

dri  dr2  d9  ,  .  .  ... 

ri— -  r2-^  =  ror2— sm0  -  rocos0  (A6) 

Solving  equations  (A5)  and  (A6)  simultaineously  for  dri/dt  and  dr2/dt  provides  the 
velocity  of  the  target  in  the  directions  given  by  ri  and  r2  respectively: 

dfi  f2  ( cwd  d9  .  „  A  , . 

= -jT  =  — T — (:7-~+r-or2— sm0-rocos0)  (A7) 

dt  ri  +  r2  \2Trfc  dt  ) 
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Figure  Al:  Geometry  of  the  velocity  of  the  target  broken  into  the  components  along  the 
lines  ri  and  r2. 

df2  —ri  ( cwd  dd  \ 

Note  that  these  equations  depend  on  dO/dt,  the  rate  of  change  of  the  look  direction,  a 
value  which  typically  requires  a  number  of  radar  “scans”  to  estimate. 

The  total  velocity  of  the  target  is  given  by  summing  these  separate  components.  This 
is  illustrated  in  figure  Al.  Using  similar  rules  of  geometry  that  were  applied  in  previous 
equations,  the  total  target  velocity,  v,  in  the  direction  is: 


V  =  +  V^2  +  2VriVr2  COs(0') 

(A9) 

and 

V  1 

0v  =  arcsin  —  sin(0') 

,  V 

(AlO) 

where  6  is  the  angle  between  ri  and  r2  and  is  calculated  using: 

9'  =  arcsin  —  sin  6 

Ln 

(All) 

A. 3  Moments  of  the  ambiguity  function  of  white  noise 

Section  4.2  gives  two  statistical  properties  of  the  ambiguity  function  when  the  signal 
of  interest  is  white  Gaussian  noise:  the  function’s  mean  and  variance.  This  section  shows 
how  these  results  were  attained. 

The  ambiguity  function  for  discrete  signals  is  defined  by  equation  (10).  Properties  of 
the  model  for  x[n]  are  given  in  section  2.2.  From  these,  it  can  be  shown  that: 

1  ^ 

=  -X!E{a:[n]a:*[n  +  r]e-*^"} 

n=l 

=  I  for  (T,d>)  =  (0,0) 

I  0  otherwise 


64 


DSTO-TR-0809 


E{|x(7-,<^)P}  =  Y,E{x[ni]x*[nx  +  T]e 


ni=l  n2=l 


(x[ni]x'[ni  +  T]e-‘*"'y) 

=  4E  EEW"lK[ni+T]e-‘*"‘} 

fll=l  712=1 

E{x*[n2]«[n2  + 

+  ]^  I]  E  E{a;[ni]x*[n2]} 

m=i  Ti2=i 

E{a;*[ni  +  r]e“*'^”^a;[n2  + 


When  (t,  0)  =  (0,0),  this  simplies  to: 


E{|x(0,0)p}  =  -i  Y.  + 


tll=l  712=1 

=  (1  +  ^)^* 


For  (r,  <!>)  pairs  away  from  the  origin: 


771=1  772=1  77=1 

=  Irr^ 

N  ^ 


E{|x(r,<^)|  }-|  ^^4^  Otherwise 


E{|x(T,</.)in  = 


,  iV  JV  JV  JV 

jfi  E  E  E  E 

ni=l  712=1  713=1  714=1 

E{a:[ni]x*[ni  +  r]e~*'^^a;[n2  +  r]e*‘^^a:*[n2] 

x[n3]a:*[n3  +  r]e“*'^^a:[n4  +  r]e*‘^^x*[n4]} 

^  N  N  N  N 

E  E  E  E 

7ll=l  772=1  713=1  7l4=l 

^E{x[ni]ar*[ni  +  T]e“*‘^^}E{a:[n2  +  T]e*‘^^x*[n2]} 
E{a:[n3]a;*[n3  +  T]e“*'^®}E{a;[n4  +  T]e*‘^^a;*[n4]} 

+  E{a:[ni]a;*[ni  +  T]e“*‘^^}E{x[n2  +  T]e*'^^x*[n2]} 
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E{a;[n3]s*[n4]}E{x[n4  -I-  T]e^H*[nz  +  r]e“**^3 
+  E{a:[ni]a;*[ni  +  r]e-*^i}E{a;[n4  + 
E{s[n3]a;*[n2]}E{a;[n2  +  T]e*'^*s*[n3  + 

+  E{a:[ni]a;*[ni  +  r]e-*^i}E{a:[n4  +  x*[n2]} 

E{a:[n3]a:*[n4]}E{x[n2  +  r]e*^=®a;*[n3  +  T]e-^'^^} 

+  E{a;[ni]x*[ni  +  T]e-*^"i}E{x[n2  +  T]e*^^a;*[n4]} 

E{x[n3]a;*[n2]}E{ar[n4  +  T]e^‘^^x*[n3  +  r]e~*'^"3j, 

+  E{a;[ni]a:*[ni  4-  r]e-*^i}E{rc[n2  + 

E{a:[n3]a:*[n3  +  T]e-^^^}E{x[n^  +  T]e*^ia;*[n2]} 

+  E{a;[ni]a:*[n2]}E{a;[n2  +  r]e*‘^^x*[ni  +  T]e'“®'^i} 
E{x[n3]rc*[n3  +  T]e-^^^}E{x[n^  +  Ty'^^x*[n^]} 

+  E{a;[ni]x*[n2]}E{a;[n2  +  TY'^^x*[ni  +  T]e“®^i} 
E{a;[n3]a;*[n4]}E{a;[n4  +  r]e®'^"ia;*[n3  +  r]e“®'^®} 

+  E{a;[ni]a:*[n2]}E{a;[n2  +  r]e*'^^x*[n3  +  r]e“*'^®} 

E{x[n3]x*[ni  +  r]e-*^i  }E{a;[n4  +  Ty^^x*[ni]} 

+  E{ar[ni]a:*[n2]}E{a:[n2  +  r]e*'^^a;*[n3  + 

E{a;[n3]a;*[n4]}E{x[n4  +  T]e*'^^a:*[ni  +  T]e~*'^^} 

+  E{a:[ni]a:*[n2]}E{a:[n4  +  r]e‘^ix*[n3  +  rje"’^®} 
E{a:[n3]x*[ni  +  T]e-*^i}E{x[n2  +  T]e*^2a:*[n4]} 

+  E{a;[ni]a:*[n2]}E{a:[n4  +  r]e*^ix*[ni  + 

E{a:[n3]®*[n3  +  T]e-^'^^}E{x[n2  +  T]e^'‘^^x*[n4\} 

+  E{a:[ni]x*[n3  +  r\e-^’>^^}E{x[n^  +  T]e*^ix*[n4]} 

E{a;[n3]a:*[n2]}E{x[n2  +  T]e^^^x*[ni  +  r]e“*'^^} 

+  E{a;[ni]x*[n3  +  T]e-*^'‘3}E{x[n4  +  T]e*^ia;*[n2]} 
E{a:[n3]a;*[n4]}E{a;[n2  +  T]e®'^"=*a;*[ni  + 

+  E{x[ni]x"[n3  +  T]e-*^®}E{x[n2  +  TY^H*[n2]} 
E{x[n3]x*[ni  +  T]e-*^"i}E{x[n4  +  r]e*^ix*[n4]} 

+  E{x[ni]x*[n3  +  T]e-*^3}E{x[n2  +  Ty'^^x*[n2]} 

E{x[n3]x*[n4]}E{x[n4  +  T]e*‘^^x*[ni  + 

+  E{x[ni]x*[n3  +  r]e“®'^®}E{x[n2  +  x*[n^} 

E{x[n3]x*[ni  +  T]e-*^i}E{x[n4  +  r]e*^ix*[n2]} 

+  E{x[ni]x*[n3  +  T]e-*^®}E{x[n2  +  T]e*^2x*[n4]} 

E{x[n3]x*[n2]}E{x[n4  +  T]e*‘^ix*[ni  + 

+  E{x[ni]x*[n4]}E{x[n2  +  T]e®‘^"^x*[ni  + 
E{x[n3]x*[n2]}E{x[n4  +  T]e*'^^x*[n3  + 

+  E{x[ni]x*[n4]}E{x[n2  +  T]e®'^^x*[ni  + 

E{x[n3]x*[n3  +  r]e-*'^"«}E{x[n4  +  T]e^'^^x*[n2]} 

+  E{x[ni]x*[n4]}E{x[n4  +  Ty‘^^x*[m  + 
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E{a:[n3]a;*[ni  +  T]e“*'^^}E{a;[n2  +  r]e*^^x*[n2]} 

+  E{x[ni]x*[n4]}E{x[n4  +  T]e*'^^x*[ni  + 
E{x[n3]x*[n3  +  r]e“*'^^®}Ex[n2  +  T]e*'^*x*[n2] 

+  E{x[ni]x*[n4]}E{x[n2  +  T]e*'^^x*[n3  +  T]e“®^®} 
E{x[n3]x*[ni  +  r]e“*'^^}E{x[n4  +  r]e*'^‘x*[n2]} 

+  E{x[ni]x*[n4]}E{x[n2  +  r]e‘‘^^x*[n3  + 
E{x[n3]x*[n2]}E{x[n4  +  r]e*'^”^x*[ni  + 

When  {T,(f>)  =  (0,0),  this  simplies  to: 

E{|X(0,0)|^}  =  ^al{N^  +  6N^  +  8N^  +  9N) 

And  for  (r,  (p)  pairs  away  from  the  origin: 


So: 


E{|x(T,.?)in  =  i^S(3JV=+6iV) 


E{|x('r,<^)|'‘} 


i^  +  W  +  W'^  for  T  =  (/)  =  0 


otherwise 
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Appendix  B:  Matlab  Code 


Two  equations  which  were  commonly  referenced  throughout  this  report  were  equations  (6) 
and  (10),  which  calculated  the  matched  filter  output  and  the  ambiguity  function  respec¬ 
tively.  This  appendbc  provides  Matlab  code  which  implements  these  equations. 

It  is  worth  noting  that  these  functions  are  very  similar.  This  is  not  surprising  given  the 
similarity  of  the  equations  and  that  both  operations  are  effectively  correlating  one  signal 
with  another  delayed  in  time  and  shifted  in  Doppler. 


B.l  Matched  filter  receiver 


function  M  =  match(xr,xt,fs,tau,phi,win) 
y#  Matched  filter  radar  receiver 

•/. 


y.  M  =  MATCH (xr,xt,fs,tau, phi, win) 

y. 


y.  xr 

y.  xt 

y. 

y.  fs 

y.  tau 

y. 

y. 

y. 

y. 

y. 

y.  phi 

y. 

y. 

y. 

y.  win 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 


Received  signal. 

Transmitted  signal.  This  must  be  the  same  length 
as  xr. 

Sampling  frequency  of  xr  and  xt  (in  Hz) . 

Vector  of  time  delays  (in  secs)  at  which  to 
evaluate  the  matched  filter. 

This  should  be  an  evenly  spaced  vector  with  the 
spacing  between  elements  being  a  multiple  of  1/fs. 
If  tau  is  the  empty  matrix,  or  is  not  specified, 
the  default  is  CO:l/fs:length(xr)/fs] . 

Vector  of  frequency  shifts  (in  Hz)  at  which  to 
evaluate  the  matched  filter. 

If  phi  is  the  empty  matrix,  or  is  not  specified, 
the  default  is  [-fs:fs/50:fs] . 

VIhat  window  function  is  used  when  processing. 

If  not  specified  or  is  the  empty  matrix,  no  window 
(effective  rectangular  or  boxcar  window)  is  used. 

If  win  is  a  vector  of  the  same  length  as  xr  and  xt, 
it  is  used  as  the  window  function, 
win  can  be  one  of  the  following  strings,  in  which 
case  the  window  function  is  calculated  and  used: 
^banning’,  ’hamming’,  ’blackman’ ,  ’boxcar’ 


y.  The  result  of  the  matched  filter,  M,  is  a  matrix  of 
It  length(tau)  rows  and  length(phi)  columns.  The  rows 
y,  correspond  to  the  possibility  of  targets  at  the  various 
It  ranges  which  (through  the  bandwidth  used  and  radar 
y,  configuration)  correspond  to  the  values  in  tau.  The 
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%  columns  correspond  to  the  radial  velocity  of  these  targets 
y.  which  correspond  to  the  frequency  shifts  phi. 

y. 

y.  Ezaiiq>le  of  use: 
y.  fs  =  1280; 

y,  xt  =  sin(2*pi*55*[0:l/fs:127/fs])  ; 

y,  xr  =  [zerosd.lO)  xt(l:118)]; 

y,  tau  =  C0:l/fs:20/fs]  ; 

y,  phi  =  [-fs/2:fs/4:fs/2]  ; 

y,  M  =  match(xr,xt,f s,tau,phi)  ; 

y,  mesh(phi,tau,M) ,  axis  tight 

y. 

y.  By  Maurice  Ringer,  1998 

error (nargchk ( 3 , 6 , nargin) ) ; 

nx  =  length (xr) ; 
if  nx‘'=length(xt) 

error (’xt  and  xr  must  be  the  same  length’) 
end 

if  nau:gin<4,  tau  =  []  ;  end 
if  nargin<5,  phi  =  [] ;  end 
if  nargin<6,  win  =  [] ;  end 

if  isempty(tau) ,  tau  =  [0:l/fs:nx/fs] ;  end 
if  isen^ty(phi) ,  phi  =  [-fs:fs/50:fs] ;  end 
if  iseii5>ty(win) ,  win  =  ones(nx,l);  end 
if  isstr(win) 

win  =  lower (win) ; 
if  strcmp (win , ’ banning ’ ) 
win  =  0.5*(1  -  cos(2*pi*[l:nx] ’/(nx+1))) ; 
elseif  strcmp (win, ’hamming’) 

win  =  0.54  -  0.46*cos(2*pi*[0:nx-l] ’/(nx-1)) ; 
elseif  strcmp (win , ’ blackman ’ ) 

win  =  (0.42  -  0.5*cos(2*pi*(0:nx-l)/(nx-l))  +  ... 
0.08*cos(4*pi*(0: nx-1) /(nx-1))) ’ ; 
elseif  strcmp (win, ’boxcar’) 
win  =  ones(nx,l); 
else 

warning  (’Funny  window  Eirgument  -  no  window  used’); 
win  =  ones(nx,l) ; 
end 
end 

ntau  =  length(tau) ; 
nphi  =  length (phi) ; 
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xr  =  xr(:) .♦win; 
xt  =  xt ( : ) . ♦win ; 
phi  =  phi(:) ; 


M  =  zeros (ntau.nphi) ; 
for  itau  =  l:ntau, 
n  =  round(tau(itau)^f s) ; 
if  (abs(n)<nx) 
if  (n>=0) 

temp  =  xr(l+n:nx)  .♦  conj(xt(l:nx-n)) ; 
t  =  [n/fs:l/fs:(nx-l)/fs]; 
else 

temp  =  xr(l:nx+n)  .♦  conj (xt(l-n:nx)) ; 
t  =  [0:l/fs: (nx+n-l)/fs] ; 
end 

M(itau,:)  =  (exp(-i^2^pi^phi^t)  ♦  temp).’; 
end 
end 

M  =  abs(M) . “2; 

M  =  M  /  maix(max(M)) ; 


B.2  Ambiguity  function  generation 


function  X  =  amb(z,t ,tau,nu) ; 
y.  Radar  Ambiguity  function 

7. 

7,  X  =  AMB(z,t,tau,phi) 

7. 

7t  z  :  Analytic  input  signal.  Values  for  z  outside  time  t  are 
7i  assummed  to  be  zero. 

7.  t  :  Timebase  for  the  input  signal  (in  seconds) .  This  must 
It  be  a  vector  of  equally  spaced  numbers . 

7.  tau  :  Vector  of  time  delays  (in  seconds)  to  evaluate  X  at 
V,  phi  :  Vector  of  Doppler  shifts  (in  Hz)  to  evaluate  X  at 

7. 

7t  X  :  Resulting  radar  ambiguity  function.  This  is  a  real 
7t  matrix  of  size  (length(tau)  .length (phi)) 

7. 

7i  Example  of  use: 

7i  t  =  [0:63]’;  z  =  exp(i^2^pi^4/64^t) ; 

7.  tau  =  [-64:4:64]’;  phi  =  [-0.16:0.016:0.16]’; 

7.  X  =  amb(z,t,tau,phi)  ; 

7.  surf  (phi,  tau,  lO^loglO(X));  axis  tight; 


70 


DSTO-TR-0809 


•/. 

7,  This  code  iir5)lemeiits  equation  4-2a  (pg  59) ,  from 
7,  "Radar  Signals:  An  Introduction  to  Theory  and  Application" 
7.  by  C  E  Cook  and  M  Bemfeld,  1967,  Academic  Press 

7. 

7t  Maurice  Ringer,  1998 

error (nargchk (4,4, nargin) ) ; 

z  =  z(:); 
t  =  t(:); 
tau  =  tau( :) ; 
phi  =  phi ( : ) ; 
t_sp  =  t(2)  -  t(l) ; 

nt  =  length(t) ; 
ntau  =  length(tau) ; 
nphi  =  length (phi) ; 

if  nt*'=length(z) 

error (’ length  of  arguments  z  and  t  must  be  the  same’); 
end 

X  =  zeros (ntau, nphi) ; 
for  itau  =  l:ntau, 

n  =  round(tau(itau)/t_sp) ; 
if  (abs(n)<nt) 
if  (n>=0) 

temp  =  z(l:nt-n)  .*  conj (z(n+l :nt)) ; 

X(itau,:)  =  temp.’  *  exp(-i*2*pi*t(l:nt-n)*phi’) ; 
else 

temp  =  z(-n+l:nt)  .♦  conj (z(l:nt+n)) ; 

X(itau,:)  =  temp.’  *  exp(-i*2*pi*t(-n+l:nt)*phi’) ; 
end 
end 
end 

X  =  abs(X)  .  '‘2; 

X  =  X  /  m8Lx(max(X)) ; 
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This  report  investigates  the  suitability  of  a  number  of  broadcast  band  transmitters  of  opportunity  for  a 
passive  radar.  Specifically,  this  report  reviews  the  design  of  active  radar  waveforms,  measures  emissions 
for  representative  transmitters  of  opportunity,  and  assesses  the  suitability  of  these  signals  for  various 
types  of  radar  observation. 

This  report  shows  that  it  is  indeed  possible  to  construct  the  suggested  passive  radar,  although  perfor¬ 
mance  is  limited  and  depends  largely  on  other  parameters,  such  as  receiver  bandwidth,  antenna  pattern 
and  observation  time. 
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